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^ 0. Preface 

U 

^ The idea for a cluster algebra and quiver package in the open-source computer 

algebra system sage was born during the sage days 20.5 which were held at the 
Fields Institute in May 2010. The purpose of this package is to provide a platform 
to work with cluster algebras in graduate courses and to further develop the theory 
by working on examples, by gathering data, and by exhibiting and testing conjec- 
tures. In this compendium, we include the relevant theory to introduce the reader 
to cluster algebras assuming no prior background; this exposition has been written 
to be accessible to an interested undergraduate. 
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Version 1.0 of the software package and this compendium is the result of many 
discussions on mathematical background and on implementation algorithms, and of 
many, many hours of coding. It is part of the sage-combinat project |'Sage Comb| . 
For more information on sage, in particular on a detailed description how to 



install the program, we refer the reader to http : / /www . sagemath . org| |Sage| ; for 
more on the sage-combinat project, see http://wiki.sagemath.org/combinat] 

Throughout this compendium, we include examples that the user can run in the 
sage-notebook or sage command line. The package provides as well an interactive 
mode for the sage-notebook as shown in Figure [T] at the end of Section |3] We 
will close with a detailed description of the data structures and methods in this 
package. We follow the usual sage convention of indexing all lists starting at zero. 

Currently, installing the sage-combinat queue is a necessary require- 
ment for working with the cluster algebra and quiver package. In order 
to install the sage-combinat queue, you have to, after installing sage, run the 
command 

# ./sage -combinat install 

on the Unix command line. Once the sage-combinat branch is created, one can 
use the command 

# ./sage -combinat update 

to update to the latest version of the sage-combinat queue, or one can use the 
command 

# ./sage -combinat upgrade 

to update to the latest version of the sage-combinat queue and to upgrade sage to 
its newest version. For more detailed explanations, please visit the sage-combinat 
wiki page. Installing the sage-combinat queue will eventually become obsolete 
after the project has gone through testing and reviewing processes, which might 
take time due to the involvedness of the algorithms (especially on mutation type 
detections) . 

This current version 1.0 should not be considered a complete, unchangeable, to- 
tally stable version. We will keep working on this project by fixing bugs, improving 
algorithms, and by adding functionalities. So it might be a good idea to update 
the sage-combinat queue once in a while, especially if you have encountered a 
problem. 

We anticipate this ongoing project being improved with feedback from users. We 
are very interested in getting any type of feedback: on ways in which the package 
has been used, on features people like or that could be done better, or requests 
for new functionalities. If you are interested in helping us make improvements or 
further develop this package, we would be happy to have you involved. 

Several other people have also worked on software for computations involving 
cluster algebras and quiver representations. Links to these are available at Fomin's 
Cluster Algebra Portal http : //www. math. Isa.umich. edu/^f omin/cluster .html. 
This software includes work of Chapoton jCha) , Dupont-Perotin |DP] , Keller |Kel| , 
and L. Williams. 
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1. Introduction 

Cluster algebras, invented by Fomin and Zelevinsky |FZ02a) . are certain commu- 
tative algebras which are isomorphic to subalgebras of the field of rational functions. 
Each cluster algebra has a distinguished set of generators called cluster variables; 
this set is a union of overlapping algebraically independent finite subsets called 
clusters, which together have the structure of a simplicial complex. The clusters 
are related to each other by binomial exchange relations. In the past ten years, such 
algebras have been found to be related to a number of other topics such as quiver 
representations, tropical geometry, canonical bases of semisimple algebraic groups, 
total positivity, generalized associahedra, Poisson geometry, and Teichmiiller the- 
ory. 

Usually, when one defines an algebra A, one describes it by writing down the 
generators and relations of A. Instead, when working with a cluster algebra, only 
a finite set of generators are provided at first, along with combinatorial data that 
allows one to algebraically construct the rest of the generators by applying a se- 
quence of exchange rules. With this definition in mind, a seed for a cluster algebra 
^ is a pair (x,B), where x denotes the initial cluster, and B denotes an exchange 
matrix (or _B-matrix|^ Here, the cluster x consists of exchangeable generators, 
known as cluster variables and non-exchangeable generators, known as coefficients 
or frozen variables. 

One of the simplest families of cluster algebras are those which are coefficient- 
free and rank two. Such algebras are parametrized by two positive integers (6, c), 
and the associated cluster algebra A{b, c) is defined to be the algebra generated by 
the set {xn\n&, where for n ^ {0, 1}, 

Xn = II n IS even, and ii n is odd. 

Xn-2 Xn-2 

These are implemented in sage, for example (letting b ~ 2, and c = 3) as 

sage: S23 = ClusterSeedC ['R2' , [2,3] ,2] ) ; S23 

A seed for a cluster algebra of rank 2 of type ['R2' , [2,3] ,2] 

Here, 'R2' refers to "rank 2", [2,3] gives the parameters. For an explanation of 



the final 2, we refer to Section 6.2 Notice that if instead we let 6 = 1 and c = 1, 
we obtain 

sage: Sll = ClusterSeedC ['R2' , [1 , 1] ,2] ) ; Sll 

A seed for a cluster algebra of rank 2 of type ['A',2] 



Technically, this is the definition for a seed of a cluster algebra of geometric type. We give a 
more general definition of cluster algebra seeds in the next section. 



4 



GREGG MUSIKER AND CHRISTIAN STUMP 



We will see more examples of this phenomenon in a moment, but the point is that 
when (6, c) = (1, 1), the associated cluster algebra is of "type A2" . 

Let us keep working with the cluster seed 511 at the moment. We can see the 
i3-matrix and initial cluster corresponding to this seed quite easily. 

sage: Sll . cluster () 



[xo,xi] 

sage: Sll .bjiatrixO 



1 

-1 



Using this data, it is possible to construct the other generators of A{1, 1) by applying 
a sequence of exchanges. We define mutation in general down below. For now, let 
us mention that if we start with the initial cluster [a;o, Xi], and mutate in the 0th 
direction, we replace the xq with X2, defined as X2 — This gives us a new 

seed, whose cluster is [2:2, xq]. 

sage: Sll .mutate (0) ; Sll . clusterO 

'xi + 1 



The exchange matrix of this new seed is simply —B. 
sage: Sll .b_matrix() 




We can continue this procedure, and now mutate in the 1st direction, letting 
X3 — replace xi. 

sage: Sll.mutate(l) ; Sll . cluster () 

Xi + 1 Xo+Xl + 1 
Xo ' XoXi 

sage: Sll .b_matrix() 



Notice that after this mutation, the exchange matrix of the obtained seed is B 
again. Consequently, we can iterate this procedure, applying the mutate command 
over and over. If we want to do this more efficiently, we can as well call mutate 
with a list of indices to apply from left to right. 

sage: Sll .mutateC [0, 1 ,0, 1] ) 

If we are not only interested in the final seed, we can instead use the procedure 
mutation_sequence. Before doing that, we reset the cluster to the initial sequence 
of variables (in the initial seed). 

sage: Sll.reset_cluster() ; 
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sage : Sll .mutation_sequence ([0,1,0,1,0] ,return_output=' matrix' ) 



1 

-1 



1 

-1 



1 

-1 



sage: Sll .reset_cluster() ; 

sage: Sll.iimtation_sequence([0,l,0,l,0] , return_output= ' var ' ) 



Xl + I Xo + Xl + 1 Xo + 1 



Xo 



XqXi 



Xl 



,Xo,Xi 



Here, the first command returns the sequence of exchange matrices obtained 
from this sequence of mutations, including the initial one. Notice, the sequence 
is exactly [B, —B, B, —B, B, —B]. The second command returns the list of cluster 
variables encountered as these exchanges occur. In the rank two case, this list is 
equivalent to [x2,X3,X4,X5,xg] corresponding to the (6, c) = (1, l)-sequence 
referred to above. 

Notice, that we have already found an interesting pattern, that is after five 
exchanges, we have arrived back essentialljj^at the same seed with which we started. 
This is therefore known as a cluster algebra of finite type and finite mutation type. 
Both of these concepts will be described in more detail below. 

For our next example, we look at the (6, c) = (2, 2) case, again a rank two cluster 
algebra. 

sage: S22 = ClusterSeedC ['R2' , [2,2] ,2] ) ; S22 

A seed for a cluster algebra of rank 2 of type ['A' , [1, 1] , 1] 

Here again, notice that this specific rank two cluster algebra is recognized. In this 
case, this is our notation for a cluster algebra of affine type Ai i. We again, run 
the procedure mutation_sequence, and obtain the following: 

sage: S22.mutation_sequence( [0,1, 0,1,0] , return_output= ' var ' ) 

xf + l xf + x^+2xf + l xf+x^ + 2x^x^+3xf+2x^+3xf + l 
xo ' x^jXi ' ^0^1 ' 

3:f +3:^+23:^3: j+33:g3;f + 43:f + 33:o+63:g3;j+63:f + 33:0 + 43: j + 1 
.,4 a > 
■^0-^1 

xl°+x^+2x^xl+3xgxj+4xlx1+5x^+4x^+9xgx1 + 12xlxj + 10x1+6x^ + 12xlx^ + 10xj+4xg + 5xl + l'\ 

^0^1 J 

Unlike the previous case, the cluster variables appear to be getting more and more 
complicated, and that pattern continues. To understand these expressions better, 
we plug in the value 1 for xq and xi. 

sage: [cv. subs(S .x(0) : 1 ,S.x(l) : 1) for cv in ms] 

[2, 5, 13, 34, 89] 

From this data, one might conjecture, and it is in fact true, that the sequence 
{x„ : XnXn-2 = x^^-i + 1 and xq ^ xi ^ 1} 
is precisely the Fibonacci numbers with even index. 

^To be precise, this seed uses matrix — B (e quivalently -B^) instead of B, but these seeds are 
the same "up to equivalence", see Remark 3.4 
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It is also clear that the cluster variables x„'s obtained by an instance of the (6, c)- 
sequence are rational functions in the indeterminates Xq and Xi. More surprisingly, 
in spite of the divisions appearing, all such x„'s are actually Laurent polynomials, 
i.e. in the ring Z[xJ^, x^^]. This is actually a special case of one of the first major 
results in the theory of cluster algebras. 

Theorem 1.1 (Laurent Phenomenon |FZ02al rFZ02bj ). Given any cluster algebra 
A, which is parameterized by a choice of exchange pattern, a choice of coefficients 
(whose group ring is given as ZPj and a choice of initial cluster {xq, xi,X2, ■ ■ ■ , Xn-i\ 
of generators, then all other generators, i. e. cluster variables, are Laurent polyno- 
mials in the ring 'LV[x^^,x^, . . . ,x^^^. 

In the same paper in which Fomin and Zelevinsky prove this Laurent phenomenon, 
they made the following positivity conjecture. 

Conjecture 1.2 (Positivity Conjecture). Given a cluster algebra A with an ar- 
bitrary exchange pattern, choice of coefficients P, and an arbitrary initial cluster 
{xq, xi, . . . , a;„_i}, then every generator of A can be written in 

Z>oP[4\a;f\...,x±li]. 

In other words, the Laurent expansions for cluster variables can be written using 
positive coefficients. 

Positivity of the coefficients is significant, as it is conjecturally related to total- 
positivity properties of dual canonical bases |FZ991 IFZOO| IZel02j . Nonetheless, this 
conjecture is still open despite nearly a decade of work by many researchers proving 
it for certain families of cluster algebras. Much of this work |CZ06[ CK08, CROSl 
WoT, 'CSOT, 'DiFKOS', ' Dup09[ [Qm) IMP071 iMSOSl IMSW091 INak09l [ProOSi iSchOSl 
[ST09, SZ04, Spe07, ZelOT] has been accomplished by exploration of examples, either 
by hand or by computer. As patterns to the Laurent polynomial expansions of 
cluster variables were noticed, the positivity conjecture and explicit formulas have 
been proven for more and more cases. This software provides further tools for such 
explorations. 

2. What is a cluster algebra? 

In this section, we give a more general and complete definition of cluster alge- 
bras, and in the next one, we describe the connection between cluster algebras and 
quivers. We say that a cluster algebra A is of rank n if ^ is subalgebra of an 
ambient field F isomorphic to a field of rational functions in n variables. Algebras 
are typically defined by generators and relations, but in the case of cluster alge- 
bras, instead of being handed all the generators at once, you are instead handed a 
distinguished set of n of them along with a constructive algorithm that can be used 
to obtain a complete set of generators. Note, that in general, a cluster algebra is 
infinitely-generated, however, any element of this distinguished generating set can 
be reached in finite time. 

This distinguished generating set is called the set of cluster variables, the first 
n of which are known as the initial cluster variables. Any set of algebraically 
independent cluster variables of maximal size has the same cardinality, namely n, 
and these n-subsets are known as clusters. Pairs of clusters x, x' whose intersection 



A COMPENDIUM ON THE CLUSTER ALGEBRA AND QUIVER PACKAGE IN Sage 7 

is of size (n — 1) are related to one another by a binomial exchange relation of the 
form 

x' = (x - {xk}) U x'k where Xkx'f, = p'^M+ + p~ M~ . 

In the second equation, and p~ belong to a coefficient semifield P (a semi- 
field is a field missing additive inverses and is subtraction-free), and Af"*", M~ are 
monomials in the elements of x — {x} which share no common factor. We let © 
denote addition in the semifield P, and multiplication is denoted in the usual way. 

Definition 2.1 (Skew-symmetrizable matrices). An n-hy-n matrix B is called 
skew-symmetrizable if there exists a diagonal integer matrix D with strictly positive 
entries on the diagonal such that DB is skew-symmetric. 

There is an algorithmic way to determine whether a matrix is skew-symmetrizable, 
and to find the diagonal matrix D, see Section [63] 



Definition 2.2 (Labeled Seed for a cluster algebra). A labeled seed for a cluster 
algebra A = .4(x, y, B) is a triple (x, y, B) where 

• X = {xq 1 } is a cluster of n algebraically independent elements 
of ambient field J- , 

• y = {j/o, VIt ■ ■ , Un-i} is an n-tuple of coefficients, elements of the semifield 
P, and 

• _B is an n-hy-n matrix that is skew-symmetrizable. 

A labeled seed can be mutated into another labeled seed {x',y',B') and all 
other clusters of A, hence all other cluster variables, can be reached by applying a 
sequence of such mutations. 

Definition 2.3 (Mutation of labeled seeds). If ^ is a cluster algebra of rank n 
and (x, y, B) is a labeled seed of A, then for any fc e {0, 1, . . . , n — 1}, there exists 
another labeled seed fik{x,y,B) = {x',y',B') = (/ifc(x), ^^(y), ^^(S)) defined as 
follows: 

The cluster x' = {xq, xi, . . . , 2?^, . . . , a;„_i} U {x'j^} where 

iVk ® l)a;fc; 

the coefficient tuple y' = {y'f^,y[, . . . ,y'^_-^) is given by 

, ^Uyr""-"\yk®ir'^' ifj^fc, 
' \i/ykiij^k 

and the matrix B' = [b'^j] is given by 

—bij if z = fc or j = fc, 

bij if bikbkj < 0, 

bij + bikbkj if bik,bkj > 0, or 
Aj - bikbkj if bik,bkj < 
We say that /ifc(x, y.B) is the mutation in the fcth direction. 

The following important observation ensures that mutation of a labeled seed is 
again a labeled seed. 
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Proposition 2.4 (Proposition 4.5 of |FZ02a| ). If B is a skew-symmetrizable ma- 
trix, then so is /ifc(-B) for < k < n — 1. 

Another helpful fact about mutation is that it is an involution, i.e. for any 
< fc < n - 1, Aifc(Aife(x, y, B)) = (x, y, B). 

Definition 2.5 (Cluster Algebras of Geometric Type). A cluster algebra of geo- 
metric type is one where the coefficient semifield P is a special choice which leads 
to a great simplification of the above formulas. In this case, 

P = Trop{uQ,Ui,U2, ■ ■ • ,Um-l) 

where the addition in Trop{uQ, ui, U2, . . . , Wm-i) is defined as 

i 3 3 

In particular, note that if yi = Y\- u^^ where all aj > 0, then ffi 1 = 1. With 
P given as this tropical semifield, the group ring ZP is simply the ring of Laurent 
polynomials Z[m^^,u^"'^, . . . ,u^_i\. 

Remark 2.6. Letting P — Trop{uo,ui, . . . ,Um^i), a labeled seed for a cluster 
algebra of geometric type is simply given as a pair (x, _B), as opposed to a triple 
(x,y,B), where x = {xq, cci, . . . , a;„_i, mq, ui, . . . , Wm-i} = {xq^xi, . . . ,Xm+n-i} 
and B is an (n + TO)-by-n matrix whose top n-by-n portion is skew-symmetrizable. 
This notation agrees with that of Section [T] 

By abuse of notation, we refer to the entire set of {n + m) variables as a cluster. 
Only the first n cluster variables are exchangeable, the last m of them are known as 
frozen variables and appear in every single cluster. These frozen variables encode 
the same data as the coefficients did in the more general case described above. The 
exchange rules for mutation instead look like the following: 

x'k^k = n ^i^" + n 

fcifc>0 bik<0 

and the mutation rule for the _B-matrix is unchanged except that we must mutate 
entires in the last m rows appropriately as well. This mutation of the last m rows 
exactly agrees with the mutation of coefficients y in the general definition. In 
particular, if we let yj = Y[a<i<m ""i*^"'^ for < j < n — 1, then we can recover the 
coefficient tuple y from the second halves of x and B. 

Remark 2.7. Since cluster algebras of geometric type are sufficient for many appli- 
cations and all of the computations currently possible in the cluster algebra package, 
we henceforth discuss the theory in terms of cluster algebras only of geometric type. 
We shall say that A = A{x, B) is a cluster algebra of rank n (with m coefficients) 
if it is a subalgebra of an ambient field F isomorphic to a field of rational functions 
in (n -|- m) variables, m of which are frozen. This is because the cluster algebra A 
is a subalgebra of ZP[xJ^, . . . ,x'^^^, and if ZP ~ Z[uq^, . . . ,u^^_i], then A can be 
thought of as a subalgebra of Z[xg ^, . . . , x^l^, wj^, . . . , u^_j], where the uf^'s are 
simply extra generators of A in addition to the set of exchangeable cluster variables. 

Note: We abuse notation and often refer to the frozen variables as coefficients; 
and we always denote frozen variables as yg through ym-i rather than the Xn+i or 
Ui notations used above. 



A COMPENDIUM ON THE CLUSTER ALGEBRA AND QUIVER PACKAGE IN Sage 9 



We close this section with some examples and more information on some basic 
commands. 

sage: B3 = matrix([[0,l,0] , [-1,0,-1] , [0,1,0]]) ; 
sage: S3 = ClusterSeed(B3) ; S3 

A seed for a cluster algebra of rank 3 

Notice that unlike the earlier examples, the description of the seed does not 
include the type. This is because the input was only the matrix, and sage will not 
attempt to recognize the type unless it is asked for by the user or by a method. 

sage: S3, cluster () 

[xo,Xx,X2\ 

sage: S3, mutate (0) 
sage: S3.b_matrix() 




sage: S3, cluster () 

"a::i + l 

We have therefore obtained a new labeled seed (x',_B') by mutating in the 0th 
direction. Note that by default, S3. mutate (0) acted on and changed the object S3 
in place. There is an option to leave S3 alone and just return the new object as a 
new output. If this behavior is desired, the command would be 

sage: S3new = S3. mutate (0,inplace=False) 

Since mutation is an involution, if we mutate again in the 0th direction, we would 
recover the original labeled seed. So we instead mutate in a different direction. 

sage: S3.mutate(l) 
sage: S3.b_matrix() 




sage: S3, cluster () 

Tl + 1 XqX2 + Xi + 1 

• ,X2 

Xo XqXi 

Let us explain why the second element (the element x'l) of this cluster is now 
3:0^:2+3:1+1 ^j^jg came from the exchange relation 

Xix'i — X2 + x'^, 

which we read off of the second column of the exchange matrix B' = ^q{B). Here 
Xq = and so we obtain the desired Laurent polynomial in terms of the initial 
cluster variables Xq, Xi, and X2 by plugging in for Xg and simplifying. 
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For one more example of the exchange relation, let us now mutate in the 0th 
direction again. This corresponds to reading the first column of B" = /ii(/io(S)) 
which gives us the exchange relation x'^ = ■ Plugging in the relevant Laurent 

polynomials for Xq and x'l, and dividing, we get a surprising cancellation and Xq is 
a Laurent polynomial: 

1\ / ^ Xi + 1\ X0X1X2 + XqX2 + Xi + 1 XqX2 + 1 



X2 



XqX2 + Xi - 
XqXi 

Using sage, we see 

sage: S3, mutate (0) 
sage : 



S3.b_matrix() 



Xq 



Xl{xi + 1) 



Xl 



sage: S3, cluster () 






-1 


1 


1 








-1 









X0X2 + 1 X0X2 +XX + 1 



, X2 



One last note about mutation, we can compress the above steps as the command 

sage: S3 = ClusterSeed(B3) ; S3 .mutate( [0 , 1 , 0] ) . 

In other words, if a list is used at the input to S3. mutate, then the seed is mutated 
to a new seed by applying the sequence of mutations in the same order as given by 
the list. 

At this point, S3 is a labeled seed with matrix B" and cluster x" as given. 
However, since a labeled seed is a choice of both an exchange matrix and a cluster, 
we also have methods to change the cluster. The first one is 

sage: S3.reset_cluster() 

This command resets the cluster to the initial cluster [xi^^xi, . . . , a;„_i] while leaving 
the exchange matrix alone. After running S3.reset_cluster(), we compute the 
exchange matrix and cluster, and obtain: 

sage: S3.b_matrix() 

i \ ll\ 

V -1 00/ 

sage: S3, cluster () 

[xo,Xi,X2\ 



A related command is S3. set_cluster() which lets the user set the initial cluster 
to be whatever they like. Note that in sage, arbitrary expressions in terms of 
indeterminates are not defined. However, integers (or even rational numbers) are 
fair to be plugged in. Additionally, if a rational function in terms of xo through Xn-i 
is desired, this can be accomplished by the commands S3 . x(0) through S3 . x(n-l) . 

sage: S3, set .cluster ( [7, 11 , 13] ) ; S3.bjiiatrix() 
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[ 1 II] 

V -1 ; 

sage: S3, cluster () 

[7,11,13] 

sage: S3.mutate([0,l,2,0]) ; SS.clusterO 

[8/11,115/77, 192/1001] 

Note that at first glance, this might seem to falsify the Laurent Phenomenon, 
but it is actually allowed because all cluster variables are supposed to be Laurent 
polynomials in terms of the initial cluster variables. Since the integers 7, 11, and 
13 are initial cluster variables, they are allowed to appear in the denominator. 

sage: SS.bjnatrixO 

»") 

V 0-10/ 

sage: S3. set_cluster ( [S3.x(0)+S3.x(l) ,S3.x(l) -2,S3.x(0)/S3.x(2)] ) 
sage: S3.cluster() 

2 Xo 

Xo + Xl,Xi, — 

sage: S3.mutate( [0,1,0,2,0]) 
sage: S3.cluster() 

xl + 1 XqxI + XoX2 + X1X2 + Xq Xoxlx2 + x\x2 + XqxI + XqX2 + XxX2 + Xq 
_Xl)+Xl' Xox\x2+x\x2 ' x1x\+Xox\ 

Again, these are Laurent polynomials in terms of the initial cluster variables ob- 
tained after setting them in this way. 

Definition 2.8 (Principal coefficients). An important cluster algebra of geometric 
type is one with principal coefficients. In this case, the initial exchange matrix B 
is 2n-by-n and where the last n rows of this matrix is a rank n identity matrix. 

Cluster algebras with principal coefficients are fundamental, because as explained 
in |FZ07) by Fomin and Zelevinsky, the formula for cluster variables in a cluster 
algebra with general coefficients (including those not of geometric type) can be 
described as a simple algebraic transformation of the formulas obtained for cluster 
variables with principal coefficients. See Theorem 3.7 of jFZ07| for more details. In 
future versions of this package, it is anticipated that working with F-polynomials 
and g-vectors will be easier with relevant methods included. 

A cluster algebra with principal coefficients can be constructed rather simply by 
the command S3.principal_extension(). Before demonstration, let us reset the 
cluster: 

sage: S3.b_matrix() 
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V -1 ; 

sage: S3.reset_cluster() ; S3.cluster() 

[10,2:1, 2:2] 

Now, we demonstrate working with principal coefficients, 
sage: SP3 = S3.principal_extension() ; SP3 

A seed for a cluster algebra of rank 3 with 3 frozen variables 
sage: SP3.b_matrix() 

/ 1 \ 
-1 1 
0-10 
1 
1 

V 01/ 

sage: SP3. cluster () 

[20, xi, 2:2, I/O, 2/1, 2/2] 

Notice, that unlike the mutate command, which is a verb, S3 is unaffected by the 
operation SP3 = SS.principal.extensionO. 

sage: S3.b_matrix() 

°°) 

V 0-10/ 

Let us try an example of mutating in this cluster algebra with principal coefficients. 
Here we use the command SP3 .mutation_sequence () instead with optional argu- 
ments return_output which is either set to 'matrix' or 'var'. 

sage: SP3.mutation_sequence( [0,1,0,2] , return_output=' matrix ' ) 
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sage: SP3.mutation_sequence( [0, 1,0,2] , return_output= ' var ' ) 

Xl + yo xpyoyi + xiX2 + 0:2^0 xpyi + X2 xoxiyoyiy2 + xpypyi + xix-z + x-zyo 

Xq ' XpXl ' Xl ' XQX1X2 

A few words about this procedure. 

(1) The command SP3 .mutation_sequence() does not affect the object SP3, 
only returns the resuhs of mutating in this order. If one wants actual seeds 
to work with rather than simply an output of matrices or cluster variables, 
one should use the option return_output='seed' (or omit this optional 
parameter since this is the default setting). 

sage: seeds3 = SP3 .mutation_sequence ( [0 , 1 , ,2] ) ; seedsS 

[A seed for a cluster algebra of rank 3 with 3 frozen variables, 
A seed for a cluster algebra of rank 3 with 3 frozen variables, 
A seed for a cluster algebra of rank 3 with 3 frozen variables, 
A seed for a cluster algebra of rank 3 with 3 frozen variables, 
A seed for a cluster algebra of rank 3 with 3 frozen variables] 

(2) With the optional parameters for returning output, the other options are 
'matrix' or 'var'. The option matrix is self-explanatory. The option 
var outputs the new cluster variable at each step. The rest of the cluster 
variables in the associated clusters are suppressed, since otherwise a lot of 
redundant information would be printed or saved. 

To return the rank (i.e. the number of exchangeable variables or columns in 
the exchange matrix B), one can simply use the command SP3.n(). The set of 
exchangeable variables can also be obtained by SP3 . exchangeable.variables () . 
To return the number of coefficients or frozen variables (also equal to the number 
of rows minus the number of columns in B), we use the command SP3.m(). The 
frozen variables are obtained by SP3 . f rozen_variables () . 

Not surprisingly, if we mutate SP3 in place with the same sequence, it equals the 
last seed returned in the sequence. 

sage: SP3. mutate ( [0, 1 ,0,2] ) ; SP3 . cluster () 



xoyi + X2 xoyoyi + xiX2 + xzyo xoxiyoyiyz + xoyoyi + xiX2 + X2yo 
, , ,yo,yi,y2 

Xl XqXi X0X1X2 

sage: SP3 == seeds3[len(seeds3)-l] 

True 



Notice that it is because of sage's indexing starting at zero that the last seed is 
indexed by len(seeds3)-l, where len stands for "length". One can also access 
the last entry by seeds3 [-1] . 

It is also rather simple to strip off the frozen variables and obtain the coefhcient- 
free cluster algebra by the command SP3 .principal_restriction() . Like the 
command S3.principal_extension(), this does not change the object in place, 
and only returns a new object where only the top half of the matrix and the first 
n cluster variables are kept. 
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Warning: Since we did non-trivial mutations before restricting, we have a la- 
beled seed where i? is a 3-by-3 matrix but the associated cluster is Laurent poly- 
nomials in terms of yo, yi, and 1/2, as well as xq, xi, and X2- 

sage: SPR3 = SP3 .principal_restriction() ; SPR3 . cluster () 

xpyi + X2 xpyoyi + X1X2 + x2yo xoXiyoyiy2 + xpypyi + xiX2 + X2yo 

Xl ' XqXi ' XpXiX2 

This causes a problem when we try to mutate: 
sage: SPR3. mutate (0) ; SPR3. cluster () 

xoyoyi + xqxi + X1X2 + x2yo xpypyi + x\X2 + x2yo 
x^yi + xpX2 ' 3:02:1 

XpXiypyiy2 + Xpypyi + X1X2 + X2yp 

XpXlX2 



since the exchange relation should have involved j/^'s but was truncated, and so we 
did not get the expected cancellation. Compare with 

sage: SP3. mutate (0) ; SP3. cluster () 

a:i-l-yo xpypy-i + X-1X2 + X2yp xpxiypyiy2 + xpypyi + X1X2 + X2yp 



Xp XpXi XpXiX2 



,yp,yi,y2 



Thus it is recommended to run SPR3.reset_cluster() or SPR3. set_cluster() 
before actually restricting to the exchangeable cluster variables by running SPR3 = 
SP3 .principal_restriction() . 



3. Using quivers as cluster algebra seeds 

In this section, we introduce a second way to input a cluster algebra seed. This 
uses the language of quivers, which is a fancy way of saying a directed (or oriented) 
graph. The term quiver originates in representation theory, where it was introduced 
by P. Gabriel at the beginning of the seventies. Gabriel wanted to emphasize the 
difference between the representation-theoretic and the graph-theoretic aspects of 
one and the same notion. For a quick introduction to quiver representations, please 
see references such as Section 5 of |Kel2] . An in-depth treatment is given, for 
example in the book by Assem, Simson, and Skowronski |ASS06| . The theory of 
quivers is an important one in representation theory, where fundamental questions 
come from studying the path algebra associated to such a directed graph. 

For our purposes, we mostly use the quivers for bookkeeping purposes and think- 
ing of them simply as directed graphs will be sufficient for most of our applications. 
In this package, a class of objects has been included as a placeholder for future 
development. For example, it is planned that in future versions of this package, 
some of the methods for quiver representations, as in preparation by Franco Saliola, 
will be available from this class as well. In the meantime, we will define what we 
need from quiver theory and describe the methods available in the current package 
as relevant to cluster algebra theory. 



A COMPENDIUM ON THE CLUSTER ALGEBRA AND QUIVER PACKAGE IN Sage 15 

Definition 3.1. A quiver Q is a directed graph. We will only work with quivers 
on a finite number of vertices and which contains no loops (1-cycles) or 2-cycles. 
However, we do allow our quivers to have multiple edges between a pair of vertices, 
but since there are no 2-cycles, this means that all edges between two vertices must 
have the same direction. In general there is no restriction against oriented cycles 
on > 3 vertices. 

Definition 3.2 (Constructing an exchange matrix from a quiver). Given a quiver 
Q on vertices vq, vi, through we let ±6ij denote the number of edges between 

Vi and Vj . We let this number be positive if the edges are oriented from Vi to Vj 
and negative otherwise. We construct Bq — [bij] as the associated n-hy-n matrix. 

Definition 3.3 (Constructing a pair-weighted quiver from an exchange matrix). 
To get a quiver Qb from an (m + n)-by-n exchange matrix B is the reverse of the 
above construction, however, there are two nuances to emphasize. 

(1) For a cluster algebra seed to correspond to a quiver, the corresponding 
matrix B must satisfy bij = —bji for all pairs < i, j < n — 1. In other 
words, the top n-hy-n portion of B must be skew-symmetric, not just skew- 
symmetrizable. Since cluster algebras for non-skew-symmetric seeds are 
also quite prevalent in the literature, our package works with a slight gen- 
eralization of quivers, which we call pair-weighted quivers. 

We do not allow parallel edges in such quivers, and instead, we label 
each directed edge as an ordered pair [bij, bji] such that the associated edge 
is oriented from Vi to Vj. Consequently, the first entry of each such pair 
is necessarily positive and the second is negative, but the direction of the 
edge must also be recorded. In the case that bji = —bij, i.e. the case of 
parallel edges, this label is simplified to be simply the positive number bij. 
We also omit the label bij = 1 when displaying graphics to make pictures 
easier to view. 

Note, that this notation differs from that in places such as |FZ03bj or 
[FSTlOj . but is necessarily for precise computations. Our notation is in- 
spired by Dlab-Ringel [DR76j and Dupont-Pcrotin jPPlOj . 

(2) lfm>0, i.e. the matrix has more rows than columns, and for any bij where 
i > n, there is no bji in the matrix and so we do not have to worry about 
checking skew-symmetry for such entries. However, such vertices Vi corre- 
spond to a frozen variable and so we designate these vertices accordingly 
as "frozen vertices" to remind the user not to mutate or apply exchanges 
at such vertices. 

Remark 3.4. This immediate connection between quivers and exchange matrices 
explains why we often consider exchange matrices up to simultaneous row and col- 
umn permutations: two quivers are considered to be isomorphic if they are isomor- 
phic as unlabeled digraphs, and this corresponds to considering exchange matrices 
up to simultaneous row and column permutations. The isomorphism reflects the 
fact that, as the cluster of an initial cluster seed {{xi, . . . , a;„}, B) is invariant under 
permuting the variable indices, the cluster algebra does not depend on the ordering 
of the vertices in the corresponding quiver. 

Definition 3.5 (Quiver Mutation). While a quiver Q can be mutated in any of 
the n directions by constructing the associated exchange matrix Bq, applying 
and then pulling back to the quiver Q^^(b) — lJ^k{Q), there is also a three step 
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process that allows for a a nice visual description of quiver mutation (in the case 
of skew-symmetric _B's). 

(1) Reverse the direction of every oriented edge incident to vertex Vk- Call the 
resulting quiver Q' . 

(2) For any 2-path — Ufe — > Vj that went through Vk in the original quiver Q, 
add a directed edge — > vj in Q' . In other words, for any pair of vertices, 
{vi,Vj}, if there are bik parallel edges from Vi to Vk and bkj parallel edges 
from Vk to Vj, then in Q', we add bikbkj directed edges between Vi and Vj. 

(3) In step 2, a 2-cycle may have been created, so the last step is to pair off 
and erase any such anti-parallel edges. 

It is an easy exercise to see that the definition of matrix mutation given 
in the previous section agrees with mutation of the quiver Qb at vertex Vk- In the 
case of a pair-weighted quiver, it is easiest to mutate the associated matrix and 
then pull-back to a pair-weighted quiver. 

Definition 3.6 (Mutation-equivalence). Two quivers Qi, Q2 are said to be mutation- 
equivalent if one can be obtained from the other by a finite sequence of mutations, 
i.e., if there exists a finite sequence ii, . . . ,ik such that /i^^. o • • • o /ij^ (Qi) — Q2- The 
collection of all quivers mutation- equivalent to a given quiver Q is called mutation 
class of Q. 

We now describe the numerous ways that a quiver can be constructed in our 
package. Firstly, a quiver can be constructed directly from an exchange matrix, or 
from a cluster seed in multiple ways. 



sage : 


B3 


= matrix ([[0,1,0] , [-1,0,-1] , [0,1,0]]) 


sage : 


S3 


= ClusterSeed(B3) 


sage : 


Ql 


= Quiver (B3) 


sage : 


Q2 


= Quiver (S3) 


sage : 


Q3 


= S3.quiver() 


sage : 


Ql 


== Q2; Q2 == Q3; Ql 

True True Quiver on i 



There are other possible constructors, such as from a directed graph: 

sage: dg = DiGraphO 

sage: dg. add_edges( [[0, 1] , [2,1]]) 

sage: Q4 = Quiver (dg) 

sage: Ql == Q4 

True 

Warning: If one uses the digraph constructor, one must follow the conventions 
for that constructor as a sage object, in particular, digraphs do not allow multiple 
edges by default. For example, to get a quiver with parallel edges, one might be 
tempted to type 

sage: dg = DiGraphO 

sage: dg.allows_multiple_edges() 

False 

sage: dg.add_edges([[0,l] , [2,1] , [2,1]]) 
sage: dg. edges (labels=False) 

[(0, 1), (2, 1)] 
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sage: Q5 = Quiver (dg) 
However, if one then asks 



sage: Ql == Q5 



True 



as the multiple copies of edge V2 — >■ vi are ignored. Instead, one should use the 
construction 

sage: dg = DiGraphO 

sage: dg. add_edges( [[0, 1, 1] , [2,1,2]]) 

sage: Q6 = Quiver (dg) 

sage: Ql == Q6; Q6 .digraphO . edges () 



Note that all quivers are actually implemented as pair-weighted quivers, i.e. as a 
labeled digraph where multiple edges correspond to a pair (6, —b) where b> 2. The 
program automatically converts the user's input with a single number indicating 
the edge label to a pair. A user can even label some edges as a single number, 
leave some edges unlabeled (as a single edge with pair- weight (1,-1)), and other 
edges as pairs; and the program will interpret this correctly. As mentioned above, 
multiple copies of an edge are ignored. More precisely, if they are given as labeled 
edges, then the label assigned is the one given to the last copy of the edge included. 

sage: dg = DiGraphO 

sage: dg. add_edges( [ [0, 1, (1,-1)] , [2,1,2]]) 

sage: Q8 = Quiver (dg) 

sage: Q6 == Q8 



A quiver can also be constructed more quickly by having sage do the intermediate 
work of constructing the digraph for you. Just simply type 

sage: Q9 = Quiver ( [ [0, 1 , 1] , [2 , 1 , 2] ] ) 

or any of the analogous constructions described above for encoding the edges of a 
digraph (although again one should include edge labels instead of multiple copies 
of edges). 

sage: Q6 == Q9 



You can also get a copy of a quiver already defined by a command such as 

sage: QIO = Quiver (Q9) 
sage: QIO == Q9 



False 



[(0,1,(2,-2)), (2,1,(1,-1))] 



True 



True 



True 



sage 
sage 



QIO. mutate (0) 
QIO == Q9 



False 



We did not emphasize it above, but a similar technique allows one to get a copy of 
a cluster seed. There is one other technique that can be used to construct a quiver. 
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or for that matter a cluster seed, but requires knowledge of quiver mutation types, 
and we leave the description of this construction to the next section. 

We now introduce some of the possible methods our package contains for working 
with quivers along with associated examples. Most importantly, to get a picture of 
a quiver, we use the show command. 

sage: Ql.showO 




sage: Ql .mutate(l) ; Ql.showO 



One quirk about the method show is that the graphic obtained can be a little 
random. If the placement of the vertices or the drawing of the graph is not optimal, 
it is recommended the user try running the show command again until the quiver 
renders in a more visually pleasing way. 

Note that just as before, the command Ql .mutateO changes the object in place. 

sage: Ql .bjnatrixO 



A nice way to visualize a sequence of quiver mutations is the use of the command 
mutation_sequence with the optional parameter show_sequence=True. Unlike the 
show command, the quivers always render in the same circular way using this 
procedure so it is easier to compare vertices to one another. 

sage: Ql.mutatioii_sequence([0,l,2,0,l,0] , show_sequence=True) 



[Quiver on 3 vertices, 
Quiver on 3 vertices, 
Quiver on 3 vertices, 
Quiver on 3 vertices, 
Quiver on 3 vertices. 
Quiver on 3 vertices. 
Quiver on 3 vertices] 





sage: Ql .mutate(l) ; Ql .bjnatrixO 
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Note that, here we are usmg the mutation.sequence on the quivers (rather than 
seeds), so the optional argument of return_output is not allowed. However, we 
can construct the associated cluster seed quite easily and then methods for viewing 
the associated quiver are still accessible, along with the other commands for cluster 
seeds. 



sage : 
sage : 



News = ClusterSeed(Ql) ; 
Q2 == Ql; NewS.showO 



Q2 = News, qui verO 



True 




sage : 



NewS .mutation_sequence( [0,1,2,0,1,0] , 

show_sequence=True , return_output= ' matrix ' ) 
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Another instructive command is DiGraph which lets the user construct the associ- 
ated labeled directed graph encoding the quiver. Since DiGraph is already a class 
of objects, this allows the user access to a variety of other methods. One can then 
reconstruct a quiver with the altered directed graph, dg whenever desired, using 
the techniques described above, i.e. Quiver (dg). 

sage: Quivs = Ql .mutation_sequeiice ( [0, 1 , 2 , 0, 1 ,0] ) 
sage: [Q.digraphO .edgesO for Q in Quivs] 

[[(0, 1, (1, -D), (2, 1, (1, -1))], 
[(1, 0, (1, -D), (2, 1, (1, -1))], 

[(0, 1, (1, -D), (1, 2, (1, -D), (2, 0, (1, -1))], 
[(0, 2, (1, -D), (2, 1, (1, -1))], 
[(2, 0, (1, -D), (2, 1, (1, -1))], 
[(1, 2, (1, -D), (2, 0, (1, -1))], 
[(0, 2, (1, -D), (1, 2, (1, -1))]] 

Thus far, the examples included have been skew-symmetric and coefficient-free. 
We close this section with some examples which require pair-weighted quivers and 
frozen vertices. 



sage : B = matrix ( [ 

(0, 1, 0, 0, 0, 1), 
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sage : 



(-1, 0, -1, 0, 0, 



(0, 1, 0, 1, 
(0, 0, -1, 



0), 
0), 
-2, 



0), 



(0, 0, 0, 1, 0, 0), 
(-2, 0, 0, 0, 0, 0)]) 
ClusterSeed(B) ; S.showO 




sage: S.mutation_sequence( [1,2,0] , show_sequeiice=True , 
r eturn_output= ' matrix ' ) 
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sage: Q = Quiver (S) 

sage: Q2 = Q.principal_extension() ; Q2; Q2.show() 

Quiver on 6 vertices with 6 frozen vertices 



^ 1 






® 







sage: Q2.mutation_sequence( [1,2,0] , show_sequence=True) 
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If we instead produce a quiver by first producing the principal_extension of the 
cluster seed, and then constructing a quiver from it, we obtain an equal quiver as 
a result. 

sage: S2 = S.principal_extension() 
sage: Q3 = Quiver(S2); Q2 == Q3 

True 

Another way to work with quivers and cluster seeds is through the interactive 
mode available through the sage-notebook. This involves a command such as 
S. interact () or Q . interact () , as shown in Figure [l] 

4. Finite type and finite mutation type classifications 

So far we have described how a cluster algebra seed can be constructed from a 
skew-symmetrizable matrix or from a quiver. The last construction that we wish 
to discuss utilizes the notion of quiver mutation types. Before, we delve more into 
the specifics of this discussion, we begin with a few theoretical preliminaries. 

Two natural questions that one can ask about a cluster algebra (or its seed) once 
the initial definitions have been given are the following: 

Given a cluster algebra A = ^(xq, So), with initial seed (xq, Bq), 

• are there a finite number of generators (cluster variables) x for A as we 
take the union of all clusters x as we mutate? 

• are there a finite number of exchange matrices B for A as we mutate into 
different seeds? 

Definition 4.1. If there are a finite number of cluster variables for A, we say that 
A is of finite type. 

Definition 4.2. If there are a finite number of exchange matrices for A, we say 
that A is of finite mutation type. 

An important theorem that greatly simplifies our notation for geometric type is 
the following theorem by Gekhtman, Shapiro, and Vainshtein: 

Theorem 4.3 (Theorem 7.4 of |GSV10j ). If A is a cluster algebra (i) of geometric 
type, or (ii) has nondegenerate exchange matrix, and (x, _B), (x',i3') are two seeds 
for A, such that cluster x' is simply the permutation a of cluster x, then the ex- 
change matrices B' and B must also he the same, up to simultaneous permutation 
of its rows and columns by the same a. In particular, the cluster determines the 
seed in the above cases. 

From this theorem, it is clear that any cluster algebra of finite type must have a 
finite number of clusters, hence a finite number of seeds and exchange matrices. 

Corollary 4.4 (Finite type implies finite mutation type). A cluster algebra of finite 
type is also of finite mutation type. 

However, the converse is false, the simplest counter-example being the rank two 
example A{2, 2) discussed in the Introduction. 

Classifying cluster algebras of finite type was one of the first natural questions 
about cluster algebras, and led Fomin and Zelevinsky to the following beautiful 
theorem. 
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S = aLEterSeed(['B',5]) 
S = 5.principa[_exterisionO 
S.iiteractQ 



Mutate at; [o]|j][2~[I][T)' 

Mutation sequence; 5^ 

Cluster variables; ^ 

B-Matrix; J 

Sliow last mutation; ^ 




Mutation sequence; 2 
Cluslier variables; 



Vj — 



Figure 1. The interactive mode of the cluster package in the sage-notebook. 
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Theorem 4.5 (Theorem 1.5 of |FZ03bj ). The following three conditions about a 
cluster algebra A — ^(xq, Bq) are equivalent: 

• Cluster algebra A is of finite type. 

• In every seed (x, J?) that is mutation- equivalent to {x.q,Bq), the exchange 
matrix B satisfies \bijbji\ < 3 for all pairs 1 < i, j < n. 

• There exists a mutation-equivalent seed (xi,Bi) such that the exchange 
matrix Bi is a skew-symmetric version of a Cartan matrix of a finite- 
dimensional Lie algebr(^ 

In particular, cluster algebras of finite type are given by the same Cartan- Killing 
classification as that describing Lie algebras via Dynkin diagrams: 

An, Bji, Cn, Dn, Eq, Ej , Eg, F4, and G2. 

Given a cluster algebra seed S for A, it therefore makes sense to ask whether 
or not S is mutation-equivalent to a seed (x, B) where the exchange matrix i? is a 
skew-symmetric version of the Cartan matrix of type An (resp. Bn, Cn, Dn, Eq, 
E-j, Eg, F4^, or 6*2)- If so, we call A a cluster algebra of mutation type An (resp. 
Bn, Cn, Dn, Eq, £'7, Eg, F4, OT 6*2). We also call all exchange matrices and the 
corresponding quivers of such a cluster algebra of mutation type An (resp. i?„, C„, 
Dn, Eq, E-j, Es, F4, or 6*2). 

Our program has algorithms for identifying mutation types of exchange matrices 
and quivers. In the cases of the exceptional types, Eq, Ey, Eg, F4 and G2, it is 
sufficient to hard-code a catalog of the mutation classes. This is done to avoid 
recomputing the mutation class whenever checking a mutation type. In classical 
types however, the parameter n can be any positive integer, and we instead utilize 
theoretical results of |CCS06] (type An), [Stull| (types i?„ and C„), and ^VatOS] 
(type Dn) to identify them for any rank n. 

Recall that a quiver (resp. pair-weighted quiver) encodes the same information 
as a skew-symmetric (resp. skew-symmetrizable) matrix. To avoid duplication of 
data types, we have introduced a new class of objects known as quiver mutation 
types. Note that these can be implemented with or without brackets. 

sage: QMl = QuiverMutationTypeC ['A' ,5] ) 
sage: QM2 = QuiverMutationTypeC ' A' , 5) ; QMl == QM2 

True 

sage: QMl 

['A', 5] 

sage: type (QMl) 

(class 'sage.alLcmdline.QuiverMutationTypeJrreducible') 
sage: QMl .bjnatrixO 

/ 1 \ 
-10-10 

1 1 

0-10-1 
\ 1 / 



Given a Cartan matrix A, we make a skew-symmetric Ba by replacing the 2's on the diagonal 
with O's, and picking a bipartite coloring of the Dynkin diagram associated to A so that bij = \aiA 
if directed edge Vi — > Vj would go from white to black, and bij = — |oij| otherwise, see Section 5 
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sage: Quiv = QMl . standard_quiver () ; Quiv 

Quiver on 5 vertices of type ['A', 5] 
sage: Quiv. show () 




sage : 



QM2 = QuiverMutationTypeCBC ,6,1) ; QM2 

['BC ,5,1] 



sage: QM2.bjnatrix() 








1 














1 \ 




-1 





-1 



















1 





1 



















-1 





-1 



















1 





2 



























v 


-2 

















/ 



sage: QM2.standard_quiver() .showO 




Each quiver mutation type has a number of attributes and methods associated to it. 
We aheady saw an example of two key methods: bjnatrix and standard_quiver, 
i.e. each quiver mutation type object encodes a specific canonical exchange matrix 
and the associated pair-weighted quiver. This characterizes only one representative 
out of the relevant possible mutation-class, but it is enough data to determine the 
appropriate cluster algebra seed up to mutation-equivalence. We hard-coded these 
representatives so that the associated quiver is an oriented Dynkin diagram such 
that each vertex is a sink or source. For future reference, such a quiver and seed is 
known as bipartite. 

More generally, each of these representative quivers are trees and acyclic. Be- 
cause of results from representation theory and otherwise, there are a number of 
results in cluster algebra theory that hold when the associated quiver is bipartite 
(resp. a tree or acyclic), but the result is incorrect, or a proof is unknown when the 
quiver lacks the relevant property. Here are some examples: 



Theorem 4.6. |Nak09j // a cluster algebra A is given by a seed that is mutation- 
equivalent to one which is skew- symmetric and bipartite, then all cluster variables 
of A have positive expansions as Laurent polynomials |^ 



^By theorems of Fan Qin [Qm| and an updated version of INakOQI , positivity has been proven 
for all skew-symmetric acyclic seeds. 
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Theorem 4.7 (Proposition 9.2 in jFZ03b| ). If Q is a quiver that is a tree as an 
undirected graph then Q is mutation-equivalent to any Q' where Q' has the same 
underlying undirected graph as Q but the edges of Q' are oriented arbitraril'^ 

Theorem 4.8 (Corollary 1.21 in |BFZ05j l. Let A = A{x,B) be a cluster algebra 
where B corresponds to an acyclic seed. Let x\ denote the unique element in cluster 
which is not contained in x. Then we have the following: 

• A is finitely generated by the set x — {xi,x[, . . . , a;„, x'„}, 

• The standard monomials (those not containing the factor Xix[ for any i £ 
{0, 1, . . . , n — 1}J in X form a ZP-basis of A, and 

• The binomial exchange relations involving Xix[ on the left-hand- sides gen- 
erate the ideal of relations among the generators x- 

Because of the importance of these properties, and other related ones, there are 
methods to check whether a given cluster seed, quiver, or quiver mutation type 
satisfies them: 

is_finite(), is_iimtation_f initeO , is_bipartite () , is_acyclic () , . . . 

There are a few other checks that we have not explained yet, but we will provide 
an annotated list of all of the checkable properties in Section |6] 

sage: QMl.propertiesO 

['A', 5] has rank 5 and the following properties: 

- irreducible : True 

- mutation finite: True 

- simply-laced: True 

- skew-symmetric : True 

- finite: True 

- affine: False 

- elliptic: False 

sage: QM2.properties() 

['BC, 6, 1] has rank 7 and the following properties: 

- irreducible : True 

- mutation finite: True 
- simply-laced: False 

- skew-symmetric: False 

- finite: False 

- af f ine : True 

- elliptic: False 

Most importantly, our program allows the user to construct a cluster seed or quiver 
by using a quiver mutation type. The associated quiver is the standard quiver that 
is hard-coded as a representative for each type; and the associated cluster seed is 
obtained from this choice of quiver. 

sage: ClusterSeedC ['A' ,5] ) 

A seed for a cluster algebra of rank 5 of type ['A', 5] 
sage: ClusterSeedC ['BC ,6,1] ) 

^Notc: this list of mutation-equivalent quivers is not exhaustive, for example a quiver of type 
is both mutation-equivalent to any orientation of a path on three vertices; or to an oriented 
triangle. 
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A seed for a cluster algebra of rank 7 of type ['BC, 6, 1] 
sage: QuiverC ['A' ,5] ) 

Quiver on 5 vertices of type ['A', 5] 
sage: Quiver ( ['BC ,6, 1] ) 

Quiver on 7 vertices of type ['BC , 6, 1] 

4.1. Finite mutation type classification. We now describe theoretical results 
regarding the classification of cluster algebras of finite mutation type. Again, we 
use the notation of pair-weighted quivers so our descriptions of some of the re- 
sults will difi^er slightly from the work of Felikson-Shapiro-Tumarkin TSTIO]. Our 
story begins however with Felikson-Shaprio-Tumarkin's first paper [FST08] which 
classified skew- symmetric cluster algebras of finite mutation type. 

Theorem 4.9 (Theorem 6.1 of |FST08j ). The following two conditions about a 
cluster algebra A — .4(xo, Bq) with skew-symmetric Bq are equivalent: 

• A is of finite mutation type, 

• A has one of the following properties: 

(1) A is of rank 2, 

(2) A is associated to a cluster algebra corresponding to a surface, or 

(3) A is one of II exceptional types Eq, Ei, Eg, affine Eq, Ey, Eg, elliptic 
E^\ E^^\ or one of two other types Xq and Xt, which were 
found by Derksen and Owen |DO08j . 

Rank two cluster algebras were already described in the introduction, and are 
clearly mutation-finite since mutation of such an exchange matrix B simply leads 
to -B. 

Describing cluster algebras of surfaces is beyond the scope of this compendium, 
however it is planned that future installments of this software will handle such 
cluster algebras and their description will be spelled out at that time. Please 
see Fomin, Shaprio, and D. Thurston's papers |FST081 IFT08| for a description or 
jMSWOQ] where Schiffler, Williams, and the first author prove positivity of Laurent 
expansions for such cluster algebras. Nonetheless, we mention here that cluster 
algebras corresponding to polygons with 0, 1, or 2 punctures, or to an annulus, can 
also be described as the skew-symmetric types A„, D„, D„, or Ar^s, respectively. 
The first two cases are of finite type and the second two are of affine type. Any other 
finite or affine type is of exceptional type or is not skew-symmetric. We illustrate 
corresponding representative quivers in the next section. 

We have met some of the eleven exceptional types before, the types Eq, E^, 
and Eg, are of finite type and thus of finite mutation type. We give representative 
quivers for the remaining eight in the next section. The affine types Eq, Ej, and Es 
each have a bipartitely oriented tree as a quiver representative; however the other 
five have no acyclic representatives. 

4.2. Skew-symmetrizable cluster algebra seeds of finite mutation type. 

In cutting edge work this summer jFSTlOj , Felikson-Shapiro-Tumarkin generalized 
their previous work to a classification including mutation-finite weighted quivers 
that are not skew-symmetric. 

Theorem 4.10 (Theorems 2.8 and 5.13 of |FST10j ). The following three conditions 
about a cluster algebra A = A{x.q,Bq) with skew-symmetrizable Bq are equivalent: 
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• A is of finite mutation type, 

• In every seed (x, _B) that is mutation- equivalent to (xqj-Bo); the exchange 
matrix B satisfies \bijbji\ < 4 for all pairs I < i, j < n. 

• A has one of the following properties: 

(1) A is of rank 2, 

(2) A is decomposable into blocks, as described in [FSTlOj . or 

(3) A is one of the 11 exceptional types in Theorem 4.9 or one of the 7 
exceptional types 6*2, -F4, -F4, V4, 1^4, ^4, md Zg. 

Remark 4.11. One can get from our notation of pair-weighted quivers to the 
notion of weighted quivers in [FSTlOj by the foUowing: if an edge of our quiver 
has the pair- weight [6, — c], then the corresponding weight in their notation is he. 
While their notation has several advantages and simplifies the statements of certain 
theorems, for computations it obscures the differences between different mutation 
classes. For example, cluster algebras of types Bn and C„ would have the same 
weighted quivers. Even though these cluster algebras give rise to the same cluster 
complexes (i.e. the clique complex induced by the graph whose vertices are seeds 
and whose edges are mutations), the Laurent expansions of cluster variables are 
quite different in these two cases. 



See Section 4.4 for 



To illustrate this example we introduce two new commands, 
details on the associated algorithms: 

1) Given a cluster algebra of finite mutation type, we can use the command 
b_matrix_class to obtain a list of all the exchange matrices that are mutation- 
equivalent to a given initial seed. To avoid extraneous duplication, we only output 
one matrix up to simultaneous permutation of rows and columns. 

For example, in the B^ versus C3 cases, notice that the list of exchange matrices 
in the respective mutation classes are negative transposes of one anothei]^ 

sage: S3 = ClusterSeedC ['B' ,3] ) ; S3.b_matrix_class() 




sage: S4 



/O 20\/0-l 

-1 1, 2 

V -1 / \ -1 1 
ClusterSeedC ['C ,3] ) ; S4.bjiiatrix_class() 




sage: S3.show(); S4.show() 



This would be clearer if we included all mutation-equivalent matrices rather than just those up 
to permutation, which could be accomplished by S3 .b_matrix_class (up_to_equivalence=False) . 
In particular the last matrices in both of these lists are negative transposes of each other if we 
also swap the first and second rows/columns. 
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(i/2| 




sage : SS.quiverO .digraphO .edgesO 

[(0, 1, (1, -D), (2, 1, (2, -1))] 
sage : S4.quiver() .digraphO . edges () 

[(0, 1, (1, -D), (2, 1, (1, -2))] 

There is an analogous command that works for chister algebras oi finite type: 

2) The command variable_class will output the list of all cluster variables 
obtained as one mutates through all mutation-equivalent seeds. 

sage: S3. variable_class() 



Xo,Xi,X2, 



Xi +1 X0X2 + 1 Xi + 1 X0X2 + Xi + 1 X0X2 + Xi 



Xo 



Xi 



X2 
V.2 I ^2 



XqXi 



X1X2 



X0X2 + Xi + 2xi + 1 X0X2 + Xi + 2xi + 1 x'l + X0X2 + 3xi + 3xi + 1 



X0X1X2 



X X"^ 



X0XIX2 



XqX% + 3x0x1X2 + xl + 2x0X2 + + 3^1 + 1 



Xi + 1 X0X2 + 1 xf + 1 X0X2 + Xi + 1 xl + XQX2 + 1 



sage: S4. variable_class() 

Xo, Xl,X2 



Xo 



Xi 



X2 



XoXi 



X1X2 



Xi + Xi+ X0X2 +2:1 + 1 X0X2 + a;i + 2xoX2 + 1 



x'fx2 



X0X1X2 

Xox\ +2:1+ X0X1X2 + x'i + 2xoX2 + a;i + 1 



Xoxlx2 

xi + X0X2 + 2x1 + 2a;oa;ia;2 + 2x1 + 2x0^2 + 2xi + 1 



2 2 

X qX -^X2 



In conclusion, even though the quivers of type B3 and C3 look quite similar and 
they have the same cluster complex, the Laurent polynomials are quite different. 
For example, the bipartite seed for a cluster algebra of type B3 leads to cluster 
variables whose numerators have degree 6, while the numerators are only of degree 
at most 4 in the case of C3. Similar phenomena happen for other dual cluster 
algebras, e.g. types i3„ versus C„ for n > 3, or pairs of seeds: (x, B) and (x, _B^). 
Here and below, we adapt the term "dual" from the notion for Kac-Moody algebras. 

Nuances like these make the non-skew-symmetric cases more difficult to analyze. 
Nonetheless, using the classification (via folding of skew-symmetric quivers) ap- 
pearing in |FST1Q| , it has been possible to include descriptions of mutation classes 
for those classes that correspond to a non-simply laced Dynkin diagram of finite or 
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afRne type, as well as the weighted quivers listed as exceptional cases in [FSTlOj . 
For the classification of non-simply laced affine Dynkin diagrams, we use the tables 
of Kac |Kac941 pgs. 53-55]. However, the notation here is not explicit enough either 
as a number of cluster algebra mutation classes are again collapsed together. We 
therefore follow notation of Dupont-Perotin |DP10j instead. The Dupont-Perotin 
notation specifies a quiver by indicating what the two ends look like, where the 
choices are that of a Dynkin diagram of type B, C or D. We say more about 
this notation in the next section. Since many users might be more familiar with 
the Kac-Moody notation, through careful coercing, if a user inputs a typical Kac- 
Moody type, it is recognized and translated into the appropriate notation that our 
software uses. 

sage: QuiverMutationType( 'C ,2) 

['B' ,2] 

sage: QuiverMutationType( 'B' ,4,1) 

['BD' ,4,1] 
sage: QuiverMutationTypeC 'C ,4,1) 

['BC ,4,1] 
sage : QuiverMutationTypeC 'A' ,2,2) 

['BC ,1,1] 
sage : QuiverMutationTypeC 'A' ,4,2) 

['BC, 2, 1] 
sage : QuiverMutationTypeC 'A' ,5,2) 

['CD', 3, 1] 
s age : Qui verMut at i onType C ' A ' , 6 , 2 ) 

['BC, 3, 1] 
sage : QuiverMutationTypeC 'A' ,7,2) 

['CD', 4, 1] 
sage: QuiverMutationTypeC 'D' ,5,1) 

['D' ,5,1] 

sage : QuiverMutationTypeC 'D' ,5,2) 

['CC ,5,1] 
sage : QuiverMutationTypeC 'D' ,4,3) 

['G' ,2,-1] 
sage: QuiverMutationTypeC 'E' ,6,1) 

['E' ,6,1] 

sage : QuiverMutationTypeC 'E' ,6,2) 

['F',4,-1] 
sage : QuiverMutationTypeC 'F' ,4, 1) 

['F' ,4,1] 

As for the finite types, our program has algorithms for identifying exchange 
matrices of affine types. In affine type A„, we have a similar coercion issue in 
the case of simply-laced affine Ar^s types where two parameters (rather than one 
parameter) is required to specify a mutation-equivalence type. This example is 
special because it is the only finite or affine type with a Dynkin diagram which is 
not a tree. Instead its Dynkin diagram is a cycle on n vertices, and here quivers Qi 
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and Q2 are only mutation-equivalent if they have the same number of edges oriented 
clockwise and the same number of edges oriented counter-clockwise. Actually, if all 
arrows are reversed, it is also the same type. The mutation classes of types A^^s 
can be classified using theoretical results in [BaslO. . 

sage: Qu = QuiverC ['A' , [2,3] , 1] ) ; Qu 

Quiver on 5 vertices of type ['A', [2, 3], 1] 
sage: Qu.showO 

sage : QuiverC [' A' , [4,1] , 1] ) . showO 
sage : QuiverC ['A' , [3,3] , 1] ) .showC) 



Notice also that the representative quiver for an affine ^r-^g-type is made as bi- 
partite as possible and that mutation type ['A' , [r,s] , 1] is coerced into type 
['A' , [s,r] , 1] when s < r. 



The remaining affine types can be found in Section 6.2 and are classified using 
resuhs in [HenQ9] and [Stull| . 

Beside the described coercions, we also include some basic coercions such as 
letting type D2 coerce into type Ai x Ai, D3 coerce into A3, C2 coerce into B2, 
small rank two examples A{b,c) coerce into A2, B2, G2, and Ai^i, and BCi for 
(&, c) — (1, 1), (1, 2), (1, 3), (2, 2), and (1,4), respectively. Here, BCi simply means 
the type ['BC ,1,1] which is a denegenerate version of the ['BC ,n,l] family of 



Dynkin diagrams used above. More technical details can be found in Section 6.2 
including other families of types and more coercions. 

4.3. Class sizes of finite and affine quiver mutation types. In this section, 
we discuss the sizes of mutation classes of finite and afhne types. Those results 
and conjectures are used to compute the size of mutation classes without explicitly 
computing the class. The class size of a cluster seed or quiver is defined to be 
the number of exchange matrices or quivers which are mutation-equivalent to the 
given cluster seed or quiver, respectively. Here, we consider seeds and quivers up 
to isomorphism. 

Theorem 4.12 (Class sizes of finite types). The number of exchange matrices or 
quivers of finite 

• type A„ [TorOSj is given by 

^2n\ / \ /2n/3 



n J \{n + l)/2j \n/3 

where the second term is omitted if {n + l)/2 is not integral and the third 
term ifn/S is not integral. 
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• type Bn or of type Cn |Stullj is given by 

1 /2n 



n + 1 \ n 

type Dn |BT09| is for n — 4 given by 6, and for n > 5, it is given by 



^ 2n \d J' 



d\n 

• types Eq, El, Eg, F^, and G2 are given by 67, 416, 1574, 15, and 2. 

Theorem 4.13 ( [BPRSlOj ). The number of exchange matrices or quivers of affine 
type Ar^s is given by 

k\r.k\s 

where (j){k) is Euler's totient function, i.e., the number of 1 < d < k coprime to k. 

Conjecture 4.14 ( [Stull| ). The number of exchange matrices or quivers of affine 

• type BBji or of type CC„ is given by 

'2n -1\ / n-1 
n-1 J ^ \n/2 - 1, 

where the second term is omitted if n is odd. 

type Dn is for n = 4 given by 9, and for n > 5, it is given by 

^2n^ 
n 



n/2j' 



2( 

where the second term is omitted if n is odd. 
type BCn is given by 

2n 

n 

• type BDn or of type CDn is given by 

J2in-l) 
\ n-l 

Theorem 4.15. The number of exchange matrices or quivers of 

• affine types Eq, Ej, Eg,, F4, and G2 are given by 132, 1080, 7560, 60, and 6. 

• elliptic types Eq^\ Ej^\ and E^^^^ are given by 49,506, and 5739. 

• the other exceptional mutation-finite types V4, W4, Xg, X7, Vg, and Zq are 
given by 7, 2, 5, 2, 90, and 35. 
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4.4. Algorithms for computing mutation classes. The four commands 

mutation_class, b_matrix_class, cluster_class, variable_class 

each utilize the auxiliary command obtained by adding _iter, which constructs an 
iterator that will run through all the objects in the corresponding mutation class. 
For quivers, there is only the method mutation_class. The first three methods are 
directly derived from mutation_class_iter, we therefore begin by describing how 
this method works. 

Note first that mutation_class_iter is, as the name already indicates, an iter- 
ator. This means that the next element is only computed if the iterator is asked to 
do so. Here is an example. One might be interested if there exists a seed or quiver 
in a given infinite mutation class having a certain property. Of course, we cannot 
test all elements, but we can construct the iterator and then let the computer run 
through the elements, constructing one after the other, and checking this property. 
If the program finds an element having the property, one could halt the process 
and return the element, together with all mutations applied to the initial element. 
If the computer keeps running, you might (or might not) get convinced that such 
an element does not exist. 

The command mutation_class_iter has five (resp. six) optional arguments if 
it is acting on a cluster seed (resp. quiver). The additional optional argument for 
quivers is data_type which is initially set to 'quiver' but can also be allowed to 
be matrix, digraph, dig6, or path. This argument does not appear in the cluster 
seed since the data type is assumed to be a cluster seed here. 

The second optional argument is depth, which is set to be 'infinity' by default 
and instructs how large a "ball" in the mutation-equivalence-class around the initial 
input is supposed to be constructed. If the cluster algebra is of finite type (resp. 
finite mutation type) however then a depth of infinity will eventually construct the 
entire mutation class, when the original input is a cluster seed (resp. quiver). 

Another optional argument is show_depth, which allows the user to print extra 
information of the actual depth, the number of constructed seeds or quivers, and the 
elapsed time. It is set to be False by default. The argument up_to_equivalence 
works differently depending on whether the input is a cluster seed or a quiver. In the 
default case True, cluster seeds are considered up to simultaneous row and column 
permutations and quivers are considered unlabeled; see Remark |3.4[ Otherwise, 
equivalence of seeds and quivers are not considered. 

sage: S = ClusterSeedC ['A' ,2] ) ; 
sage: S. cluster _class() 



Fo,a;iJ , 



a:o -I- 1 



Xl 



xo 



, Xl 



xo + xi + 1 a;o -I- 1 



XqXi 



Xl 



a;i -I- 1 a:o -I- a:i -I- 1 



xq 



XoXi 



sage : S . cluster _class (up_to_equivalence=False) 



[xo,xi] , 



Xo, 



Xl 



a;o + a;i + 1 a;i -f 1 
xoxi ' xo 



Xl + 1 



Xo 



, Xl 



Xo + 1 Xo + Xl + 1 
Xl ' XqXi 



Xl + 1 Xo + Xl -I- 1 

Xo ' xoxi 

Xl + 1 



Xl, 



Xo 



a::o + Xl + 1 Xo + 1 


XqXi 


■) 

Xl 


'xo + 1 




,xo 


, [xi,xo] 


Xl 
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The argument sink_source is set to be False by default, but if set to True, tlicn 
only mutations at sinks and sources are performed. This option is helpful for 
working with bipartite seeds or studying the BGP reflection functors on quiver 
representations. 

Finally, the last argument return_paths, again False by default, will keep track 
of the shortest mutation sequence that can be used to produce a given seed (or 
quiver) from the initial one. This data can be accessed by other commands and then 
utilized for future work. Note that such a sequence is not unique so accessing this 
shortest sequence during different computational sessions might not give the same 
result but for most purposes a single example of the mutation sequence between 
two seeds is sufficient data. 

With this iterator, one can then call mutation_class which will output the as- 
sociated list of seeds or quivers in the mutation class. However, since this output 
cannot be infinite, the argument depth cannot be infinity unless the input is of 
finite (resp. finite mutation type). The data associated to the optional arguments 
is also returned at this time. The commands b_niatrix_class and cluster_class, 
which each can only be performed on a cluster seed, work analogously. The algo- 
rithm for variable_class, which again only works on a cluster seed, requires a 
little more explanation. 

The procedure for variable_class_iter starts by running through an iterator 
for the mutation class and by yielding all found cluster variables. However, since 
the set of cluster variables is dwarfed by the number of clusters, this search-based 
algorithm is quite slow. 

On the other hand, if we are in the lucky situation that the initial cluster is 
bipartite, then we can use [FZ07t Theorem 8.8] to efficiently compute the variable 
class. 

Theorem 4.16 (Theorem 8.8 of jFZ07| ). Suppose that an exchange matrix B is 
bipartite, and its Cartan counterpart A = A{B) is indecomposable. 



1) If A is of finite type, then the corresponding bipartite belt (see Definition J^.ll) 
has the following periodicity property: the labeled seeds Y,m andT,m+2(h+2) equal 
to each other for all m £ Z,. Here, h is the Coxeter number of the corresponding 
Cartan matrix A. 

2) If A is of infinite type, then all of the elements Xi-2m, denoting the n cluster 
variables o/S2m = {{xi-2m, X2:2m, ■ ■ ■ , Xn:2m}, B) as TO ranges over the integers are 
distinct Laurent polynomials in the initial data. 

Note that in this theorem, the Cartan counterpart of B (see Section [s]) is the 
(generalized) Cartan matrix A = A{B) = (aij) defined by 




Definition 4.17. We use Sq = (^Oj^) to denote an initial bipartite seed and let 
fi+ (resp. fi-) denote the concatenation of all mutations at sources (sinks) of the 
quiver Q{Bf\ Observe that /i+(B) = ^-(B) = -B. 



Since sources and sinks are not adjacent, the factors of ^+ (resp. fj.-) commute with one 
another, hence why /^_|_ and /i_ are well-defined. 
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Define the associated bipartite belt to be the seeds S„i — (xm, (— for 
rn € Z, defined recursively by 



if r is odd 
if r is even. 



As a consequence, given an initial bipartite seed (x, B), it is sufficient to mutate 
all vertices labeling sinks in Q{B) followed by mutating all vertices labeling sources 
in Q{B), and iterate. We will get no repeats in this list and thus the most efficient 
way to obtain all cluster variables in the case of a finite type cluster algebrej^ 

Our algorithm thus first checks if the initial seed is bipartite for this reason. If 
not, it proceeds as above trying to mutate in all directions. 

It is a difficult computational problem to find a mutation sequence, if one exists, 
from an initial non-bipartite seed to a bipartite one, so it is not computationally 
feasible to use the shortcut if we do not have a bipartite seed at hand. However, 
since our proceeding is doing a search through all seeds mutation-equivalent to the 
initial one anyway as its default behavior if we get lucky and find a bipartite seed, 
the program can record this path and take advantage of this find. 

In the case that the search algorithm finds a bipartite seed, the algorithm then 
does the following procedure instead: 

1) Starts over at the initial seed. 

2) Mutates along the recorded path to get to the bipartite seed Sq- 

3) Mutate along the bipartite belt the appropriate distance from there in both 
directions (i.e. applying /i+ first or fi^ first). 

In step (3) the appropriate distance is either the period 2{h + 2) in the case of 
a cluster algebra of finite type or the depth chosen beforehand by the user. Note 
well that the meaning of depth is actually different here, as the algorithm will no 
longer spread out in all directions. Instead, the argument depth now instructs 
the computer how many iterations of the bipartite belt to use. The program will 
actually output the cluster variables found on the way to the bipartite seed Sq, as 
well as all cluster variables in the seeds {S„i : m e Z, |to| < depth}. 

Since in the case of infinite type, not all cluster variables can be reached by using 
the bipartite belt, for example even cluster variables lying in clusters two mutations 
away from the bipartite seed might not be reachable (see the bipartite ^2^2 example 
below), the optional argument ignore_bipartite_belt=False is included. If set 
to be True, the original (albeit slower) algorithm of mutating in all directions out 
to a certain depth is utilized even if a bipartite seed is found. 

sage: S = ClusterSeed(['A' , [2,2] ,1]) ; S .bjiatrixO ; S . is_bipartite() 

/ -1 -1 \ 

10 10 

0-10-1 ^^^^ 
\ 1 1 / 
sage: S.variable_class(depth=l) 

Found a bipartite seed - 

constructing the variable class into its bipartite belt. 



^ If the cluster algebra is of infinite type, one can also mutate along the bipartite belt to 
efficiently generate a large list of cluster variables but not all cluster variables are reachable in 
this way. 
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X1X3 + 1 X0X2 + 1 X1X3 + 1 X0X2 + 1 

^0; -^1; -^35 ; 5 5 ; 

Xo Xi X2 X3 

xlx'i + X0X2 + 2X1X3 + 1 X0X2 + 2X0X2 + X1X3 + 1 

X0XIX2 ' X0X1X3 

xlxl + X0X2 + 2X1X3 + 1 X^xl + 2X0X2 + X1X3 + 1 

X0X2X3 ' X1X2X3 

If we look at the output from S . variable_class(depth=2) or higher depth, we 
will see that the denominators grow larger and larger but no denominator of xqXi 
appears. Compare this output with the examples below. 

sage: S. mutate ( [0, 1] ) ; S. cluster () 

X1X3 + 1 X0X2 + X-i_X3 + 1 

5 ; X2 5 X3 

Xo XoXi 

sage: S.variable_class(depth=2, ignore_bipartite_belt=True) 

X1X3 + 1 X0X2 + 1 X1X3 + 1 

Xo,Xi,X2,X3, , , , 

Xo Xi X2 

X0X2 + X1X3 + 1 X0X2 + X1X3 + 1 X0X2 + X1X3 + 1 

XqXi ' X0X3 ' X1X2 

xf X3 + X0X2 + 2X1X3 + 1 XqxI + 2xoX2 + X1X3 + 1 

X0X1X2 ' X0X1X3 

x\x3 + X0X2 + 2X0X1X2X3 + 3xfx§ + 2xoX2 + 3xiX3 + 1 

X0X1X2X3 

5. ASSOCIAHEDRA AND THE CLUSTER COMPLEX 

Before looking at associahedra, the cluster complex and their implementations, 
we need to start with some basic background on root systems for (generalized) 
Cartan matrices. For further details, we refer to jHum72|. iKacP] . 

Definition 5.1 (Generalized Cartan matrix). An n x n-matrix A — (aij) with 
integer entries is called a generalized Cartan matrix if 

• flii = 2, 

• Qij < for j ^ j, 

• ^4 is symmetrizable, i.e., there exists a diagonal matrix D with positive 
entries such that DA is symmetric. 

A generalized Cartan matrix is called of finite type if DA is positive definite, and 
of affine type if DA is positive semi-definite. 

Recalling the definition of £?-matrices, we see that we can associate a generalized 
Cartan matrix to every _B-matrix (see |FZ03bl (1.6)]). The terms finite and affine 
come from their connections to finite and affine Lie algebras. Indecomposable gen- 
eralized Cartan matrices of finite type (resp. of affine type) classify Lie algebras of 
finite type (resp. of affine type). 

A realization of a Cartan matrix A (of finite type) is a (rational, real, or complex) 
vector space V with distinguished basis A = {at : < i < n}, and with dual space 
V* with distinguished basis = {a^ : < z < n}, together with the pairing 
(a>',aj) = a,j. For jS € V (resp. e V"^), we write [^,a,] (resp. [/?^,a^]) for the 
coefficient of in (3 (resp. in /3^). 
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Define a reflection on V by 

and define moreover, tlie Weyl group by W — {s^ : < i < n) < Aut(y) and the 
root system by 

$ = {Lj{a) : uj e W, a e A} . 

It can be shown that <f> can be written as U where 

$+ = G $ : a] > for ah a G A}, 

and $~ — {—(3 : /3 G . The elements in $ are called roots, the elements in 
are called positive roots, and the elements in A are called simple roots. 

Theorem 5.2 (Theorem 1.9 of [FZ03b| ). Let A be a Cluster algebra of finite type 
and let <l'>_i = U —A be the set of almost positive roots of the root system 
of the associated Cartan type given by the positive roots together with the simple 
negative roots. There exists a unique bisection between almost positive roots and the 
cluster variables for A for which the simple negative root ~a is mapped to Xa and, 
for positive roots, 

\ ^ Pn 



^ A 1 1 



with Pa having nonzero constant term. Here, Xa^ stands for Xi for an appropriate 
ordering A = {ag, . . . , a„_i}. 

This connection in the finite types can be used in the cluster algebra package as 
follows: 

sage: for f in ClusterSeedC ['A' ,2] ) . variable_class() : 
. . . . : print f , f . almost_positive_root () 



xo — ai 
Xl — a2 
{xi + l)/xo ai 
(a;o + l)/3;i 02 
(a;o + a;i + l)/(a::oa;i) ai + a2 

sage: f 

(a;o + a;i + l)/(a::o2;i) 

sage: root = f . almost_positive_root () ; root 

Ql + 0-2 

sage: root, parent () 

Root lattice of the Root system of type ['A' , 2] 



5.1. Generalized associahedra. In this section, we will define generalized asso- 
ciahedra and describe how they can be realized as polytopal complexes in finite 
types. We will see then how these polytopal complexes are implemented in sage. 
Generalized associahedra beyond finite type are not yet feasible as the needed tools 
to deal with infinite types are not yet developed. We start with the definition of 
generalized associahedra (not necessarily of finite type). 



A COMPENDIUM ON THE CLUSTER ALGEBRA AND QUIVER PACKAGE IN Sage 37 



Definition 5.3 (Generalized associahedron). The generalized associahedron asso- 
ciated to a cluster algebra A can be defined as the exchange graph of the mutation 
class of all cluster seeds for A. This is the unoriented graph with vertices given 
by the set of all cluster seeds, and with an edge joining two clusters if they can be 
obtained from each other by a mutation. 

Generalized associahedra reduce in classical types to known constructions, see 
e.g. |FZ03b[ Section 12]. By |FZ03b[ Theorem 1.12], a cluster seed of finite type 
is uniquely determined by its cluster, and two seeds are obtained from each other 
by a mutation if and only their clusters differ by exactly one cluster variable, see 
Theorem |4.8[ In finite types, there exist realizations as polytopal complexes, see 
|CFZQ2| . Let 5+, S'„ be the bipartition of the simple reflections S = {sa ■ ct G A} 
corresponding to the simple roots in A. This means that and S- are chosen in 
such a way that the reflections in each pairwise commute. Observe that the fact 
that all quivers of finite type are bipartite ensures that such bipartitions always 
exist. Define two piecewise linear operators and r_ on V by 




/3 \i ji — ~a for Sa G S-^ 

rises ^^l^) otherwise. 



and let 

^e*+ 

In [CFZ021 Theorem 1.1], it is shown that every (t_|_, r_) -orbit in $>_i intersects 
—A. Moreover, ai,aj G —A lie in the same orbit if and only if = —ujo{aj) 
where luq is the (unique) longest element in W. Thus, the coefficients [p^,a^] and 
[p^jOj] coincide; for (5 E $>-i, set C/j to be this coefficient. After identifying tp 
with the n-tuple {{ip, ai))o<i<n, define the half-space 

H+{P) :={^eM": <c4 

to obtain the polytopal realization of the generalized associahedron by 

Ass($) = Pi H+{(3) C M". 
I3e<s>+ 

The operators t+ and t_ are implemented in sage as operators for the root space. 

sage: S = RootSystemC ['A' ,2] ) .root_spaceO 

sage: tau_plus, taujninus = S . tau_plus_minus O 

sage: for beta in S . almost_positive_roots () : 

....: print beta, tau_plus (beta) , taujninus (beta) 

. . . . : print 



— ai, ai, —ai 
Oi, — ai, ai + Q2 

ai -I- Q!2, Q2, cti 

— a2, — a2, Q2 

Q2, Ql + Q2, —02 



sage : 



AssoA2 = Associahedron( ['A' ,2] ) ; AssoA2 
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The generalized associaiiedron of type ['A', 2] 

having 2 dimensions and 5 vertices 

sage: AssoB2 = AssociahedronC ['B' ,2] ) ; AssoB2 

The generalized associahedron of type ['B', 2] 

having 2 dimensions and 6 vertices 

sage: AssoC2 = AssociahedronC ['C ,2] ) ; AssoC2 

The generalized associahedron of type ['C, 2] 

having 2 dimensions and 6 vertices 

sage: AssoG2 = AssociahedronC ['G' ,2] ) ; AssoG2 

The generalized associahedron of type ['G', 2] 

having 2 dimensions and 8 vertices 

sage: AssoA2.showC) ; AssoB2 . showC) ; AssoC2 . showC) ; AssoG2 . showC) 




sage: AssoAS = AssociahedronC ['A' ,3] ) ; AssoA3 

The generalized associahedron of type ['A', 3] 

having 3 dimensions and 14 vertices 

sage: AssoB3 = AssociahedronC ['B' ,3] ) ; AssoB3 

The generalized associahedron of type ['B', 3] 

having 3 dimensions and 20 vertices 

sage: AssoA3.showC) ; AssoB3 . showC) 
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The associahedron of type A3 has 14 vertices (13 of which are visible, the 14th 
is the origin, which corresponds to the cluster {— ai, — a2, — 0:3}). As well, the 9 
facets corresponds to the almost positive roots, where the hyperplane Xi = c^^i 
correspond to the simple negative root — a^. Every vertex corresponds to exactly 3 
hyperplanes, and in type By,, we have 20 vertices and 12 facets, as desired. 

5.2. The cluster complex. As with associahedra, we will define the cluster com- 
plex in general and then discuss the implementation for finite types. 

Definition 5.4 (Cluster complex). The cluster complex associated to a cluster 
algebra A can be defined to be the simplicial complex with vertices being the 
cluster variables for A and with facets being the clusters. 

As we have seen, cluster variables in finite types are in bijection with almost 
positive roots. We use this description in the implementation of the cluster complex. 

sage: ClusterComplexC [' A' ,2] ) 

Simplicial complex with 5 vertices and 5 facets 
sage: ClusterComplexC ['A' ,3] ) 

Simplicial complex with 9 vertices and 14 facets 
sage: Delta = ClusterComplexC ['B' ,3] ) ; Delta 

Simplicial complex with 12 vertices and 20 facets 

In the following example, we see how we can use other sage packages to further 
study objects we work with. As the cluster complex is a simplicial complex, there 
now exists various possible methods. For example, we can compute its homology, 

sage: Delta. homologyC) 

{0 : 0,1 : 0,2 : Z} 

This is as expected, as this simplicial complex is the boundary complex of a 
triangulated polytope, and thus shellable and Cohen-Macaulay. 

6. Methods and attributes 

In this section, we describe the different classes defined in this package, and list 
their attributes and methods. For the "key" methods, we also give descriptions of 
the algorithms. 
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In general, attribute names start with an underscore to emphasize that they 
should not be used directly but only through appropriate methods. As an example, 
a cluster seed has an attribute _M in which its exchange matrix is stored and a 
method b_matrix which is used to get the exchange matrix. The difference is that 
the method returns a copy of its exchange matrix, so it is safe to work with this 
matrix and to modify it without accidentally modifying the seed itself. 

sage: S = ClusterSeedC ['A' ,3] ) ; 
sage: Ml = S._M; M2 = S.b_matrix() ; 
sage : Ml == M2 

True 

sage : Ml is M2 

False 

6.1. Skew-symmetrizable matrices. We briefly want to describe the algorithm 
used to determine whether a square matrix B is skew-symmetrizable, which also 
determines the associated diagonal matrix D in the affirmative case. It was written 
by F. Block, F. Saliola, and C. Stump during the sage days 20.5 at the Fields 
Institute, Toronto, Canada, in May 2010. 

Algorithm 6.1. Let B — {bij)i<ij<:n be the input square matrix of dimension n, 
and let D — (di)i<i<„ be the diagonal matrix with positive coefficients we want to 
construct. We use the equivalent description of skew-symmetrizablility given by the 
property 

= -djbji for all i,j. 

(1) Check if bii — for all i. If this is not the case, return False, 

(2) let k be the smallest integer such that dk is not yet determined, 

(3) set dk — 1, 

(4) for i g {1, . . . , n} such that bik 7^ and di is not yet determined, do 

(a) set di = -dkbki/bik, 

(b) if di < return False. 

(c) if any {dibij ^ ~djbji) for j such that dj is already determined, return 
False. 

(5) repeat step (4) with k given by all integers for which di was set since we 
passed step (3) the last time, 

(6) if D is not yet completely determined, goto step (2), 

(7) return D. 

6.2. QuiverMutationType. For coding reasons, we distinguish between the classes 
QuiverMutationType_lrreducible and QuiverMutationType_Reducible, but we 
refer here to both as QuiverMutationType. Objects of those types are unique, i.e., 
there exists only one object of a given quiver mutation type. 

sage: mut_typel = QuiverMutationType( ' A' , 3) 
sage: mut_type2 = QuiverMutationType( ' A' , 3) 
sage: mut_typel is mut_type2 

True 

All the data for quiver mutation types is hard-coded. In particular, this concerns 
the graphs and digraphs, and the class size. 
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To construct a quiver mutation type, the function QuiverMutationType is called. 
An irreducible quiver mutation type takes 3 parameters, the letter, the rank or 
bi_raiik, and the twist, see the description below. Those calls are best explained 
in examples. Observe that the call arguments can be also wrapped into a list or 
tuple. We suppress the output whenever the output coincide with the input. 

• finite types 

sage: QuiverMutationType (' A' ,1) ; 

sage: QuiverMutationType ('A' ,5) ; 

sage: QuiverMutationType ('B' ,2) ; 

sage: QuiverMutationType ('B' ,5) ; 

sage: QuiverMutationType ('C ,2) 

['B', 2] 

sage: QuiverMutationType ('C ,5) ; 
s age : Qui verMut at i onType ( ' D ' , 2 ) 

[ ['A>, 1], ['A', 1] ] 

sage: QuiverMutationTypeCD' ,3) 

['A', 3] 

QuiverMutationTypeCD' ,4) ; 
QuiverMutationType ( 'E' ,6) ; 
QuiverMutationType ( 'E' ,7) ; 



sage : 
sage : 
sage : 
sage : 
sage : 
sage : 

affine types 



QuiverMutationType ( 'E' ,8) 
QuiverMutationType ( 'F' ,4) 
QuiverMutationType ( 'G' ,2) 



sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

['A 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

['A 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 

sage: QuiverMutationType ( 



1] 



1] 



(1,1), 1); 
(2,4),i); 
,1,1) 
[1, 1] 

2,1) 
4,1) 

1,1) 
[1, 1] 

,2,1) 
,4,1) 
,1,1) 
,5,1) 
,3,1) 
,5,1) 
,3,1) 
,5,1) 
4,1) 
6,1) 
6,1) 
7,1) 
8,1) 
4,1) 
4,-1) 
2,1); 
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sage : QuiverMutationTypeCG' ,2,-1) ; 

hyperbolic types 

sage : QuiverMutationTypeCE' ,6, [1,1] ) ; 
sage : QuiverMutationTypeCE' ,7, [1,1]); 
sage : QuiverMutationTypeCE' ,8, [1,1]); 

mutation-finite types 



rank 2 








sage : 


Qui ver Mut at i onType 


( 


R2' , (1,1) ,2) 








['A', 2] 


sage : 


Qui ver Mut at i onType 


( 


R2' , (1,2) ,2) 








['B', 2] 


sage : 


Qui ver Mut at i onType 


( 


R2' , (1,3) ,2) 








['G', 2] 


sage : 


Qui ver Mut at i onType 


( 


R2' , (1,4) ,2) 






[ 


'BC, 1, 1] 


sage : 


Qui ver Mut at i onType 


( 


R2' , (1,5) ,2) ; 


sage : 


Qui ver Mut at i onType 


( 


R2' , (2,2) ,2) 




[' 


A 


', [1, 1], 1] 


sage : 


Qui ver Mut at i onType 


( 


R2' , (3,5) ,2) ; 


exceptional types 






sage : 


Qui verMut at i onType 


( 


V ,4,2) ; 


sage : 


Qui ver Mut at i onType 


( 


W',4,2); 


sage : 


Qui ver Mut at i onType 


( 


W ,4,-2) ; 


sage : 


Qui ver Mut at i onType 


( 


X' ,6,2) ; 


sage : 


Qui ver Mut at i onType 


( 


X',7,2): 


sage : 


Qui ver Mut at i onType 


( 


Y' ,6,2) ; 


sage : 


Qui ver Mut at i onType 


( 


Z' ,6,2) ; 


sage : 


Qui ver Mut at i onType 


( 


Z' ,6,-2) ; 



mutation-infinite types 

— infinite type E 

sage: QuiverMutationTypeCE' ,9,3) 

['E', 8, 1] 

sage: QuiverMutationTypeCE' ,10,3) ; 

sage: QuiverMutationTypeCE' ,12,3) ; 

sage: QuiverMutationTypeCAE' , (1,1) ,3) ; 

sage: QuiverMutationTypeCAE' , (1,4) ,3) ; 

sage : QuiverMutationType( 'BE' ,5,3) ; 

sage : QuiverMutationType( 'CE' ,5,3) ; 

sage : QuiverMutationTypeCDE' ,6,3) ; 

— Grassmannian types - the second parameter (a, b) must satisfy f < 
a < b and one obtains a grid graph of width a — 1 and height b ~ a~ I 

sage : QuiverMutationTypeCGR' , (2,4) ,3) 

['A', 1] 

sage : QuiverMutationTypeCGR' , (2,6) ,3) 

['A', 3] 

sage : QuiverMutationTypeCGR' , (3,6) ,3) 

['D', 4] 
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sage: QuiverMutationTypeC 
sage: QuiverMutationTypeC 
sage: QuiverMutationTypeC 

[': 

sage: QuiverMutationTypeC 
triangular types - the second 
sage: QuiverMutationTypeC 

sage: QuiverMutationTypeC 



GR' , C3,7) ,3) 
['E', 6] 
GR' , C3,8) ,3) 
['E', 8] 
GR' , C3,9) ,3) 
['E', 8, [1,1]] 
GR' , C3,10) ,3) ; 

parameter gives the size of the graph 

TR' ,2,3) 
['A', 3] 
TR' ,3,3) 
['D', 6] 
TR' ,4,3) 
', 8, [1, 1]] 
TR' ,5,3) ; 



sage: QuiverMutationTypeC 

['E 

sage : QuiverMutationType C 

type T - the second parameter gives the lengths of the three legs 



sage : 


Qui verMut at i onType 


C 


T' 


,C1,1,1),3) 








[ 


A', 1] 


sage : 


Qui ver Mut at i onType 


C 


T' 


,C1,1,4),3) 








[ 


A' , 4] 


sage : 


Qui ver Mut at i onType 


c 


T' 


, Cl,4,4) ,3) 








[ 


A', 7] 


sage : 


Qui ver Mut at i onType 


c 


T' 


, C2,2,2) ,3) 








['D', 4] 


sage : 


Qui ver Mut at i onType 


c 


T' 


,C2,2,4),3) 








['D', 6] 


sage : 


Qui ver Mut at i onType 


c 


T' 


, C2,3,3) ,3) 








[ 


E' , 6] 


sage : 


Qui ver Mut at i onType 


c 


T' 


,C2,3,4),3) 








[ 


E', 7] 


sage : 


Qui ver Mut at i onType 


c 


T' 


, C2,3,5) ,3) 








[ 


E', 8] 


sage : 


Qui ver Mut at i onType 


c 


T' 


,C2,3,6),3) 








:'E 


8, 1] 


sage : 


Qui ver Mut at i onType 


c 


T' 


,C2,3,7),3) 






[ 


'E 


, 10, 3] 


sage : 


Qui ver Mut at i onType 


c 


T' 


,C3,3,3),3) 








:'E 


6, 1] 


sage : 


Qui ver Mut at i onType 


c 


T' 


, C3,3,4) ,3) ; 


• reducible types 









sage : QuiverMutationTypeC ['A' ,3] , ['B' ,4] ) 
[ ['A' , 3] , ['B' , 4] ] 



As described in Section 4.2 one can use also Kac's classification types |Kac94| . 

Remark 6.2. Most of the above types have already been explained as Dynkin di- 
agrams, appear in Kac's list, or in the classification work of Derksen-Owen |DO08| . 
and Felikson-Shapiro-Tumarkin [FSTOS) IFSTlOj . The exceptions to these are the 
triangular seeds, Grassmannian seeds, and the "T" seeds. The first two of these 
describe a certain family of quivers that have certain shapes (as triangles and grids. 
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respectively) and correspond to certain coordinate rings of geometric objects. (See 
Examples 4.4 and 4.6 of [Kel2] or the source papers |BFZ05| and |Sco06| .) The "T" 
family are those which correspond to "Dynkin diagrams" of the shape of a T with 
a certain number of vertices on each arm and one central vertex. 

sage: ClusterSeedC ['TR' ,5,3]) .showO 
sage: ClusterSeedC ['GR' , [5,11] ,3]) .showO 
sage: ClusterSeedC ['T' , [4,4,5] ,3]) .showO 




We also illustrate a self-dual and two dual non-simply laced exceptional mutation- 
finite cases here too. 



ClusterSeedC ['X' ,6,2] ) .showC) 

S = ClusterSeedC ['W, 4, 2]); S.showC) 

S = ClusterSeedC ['¥', 4, -2]); S.showC) 



sage : 
sage : 
sage : 
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The attributes of QuiverMutationType are given by 

• .letter 

The string containing the letter (s) of the classification type. 

• _rank 

The number of vertices in the standard quiver. 

• _bi_raiik 

Is None except for afHne type A, where it denotes [a, b] with a + b being the 
rank and a < b are the number of edges in the acyclic orientation of the 
standard quiver. 

• _twist 

Depends on the type of the classification type, and can be one of the fol- 
lowing: 

— None for finite types, 

— 1 for afhne types, 

— [1,1] for elliptic types, 

— 2 for finite mutation types which are not finite or elliptic, 

— 3 for infinite mutation types. 

• -graph 

Graph representing the underlying graph of the standard quiver. 

• .digraph 

Digraph representing the underlying graph of the standard quiver. 

• -description 

The string representation of the mutation class. 

• _info 

Dictionary containing the keys 

— irreducible, 

— finite, 

— affine. 
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— elliptic, 

— simply _laced, 

— mutation_finite, and 

— irreducible_components. 

The values are True or False, except for irreducible_components which 
is a list containing the irreducible components. 

The methods of QuiverMutationType are given by 

• __eq__(self , other) 

Returns True, iff self and other represent the same quiver mutation type. 
As quiver mutation types are unique (i.e., there exists at most one object 
representing a given quiver mutation type), this method simply returns 
self is other. 

• _repr_(self) 

Returns the string representation of self. 

• plot (self , circular=False , directed=True) 

Returns a random or circular, directed or undirected plot of self. 

• show(self, circular=False , directed=True) 

Shows the plot of self. 

• rajik(self) 

Returns the rank (i.e., the number of vertices) of self. 

• coxeter_diagram(self ) 

Returns the Coxeter diagram of self 

sage: QuiverMutationTypeC ['A' ,5] ) . coxeter_diagram() 

Coxeter diagram of rank 5 
sage: QuiverMutationTypeC ['A' ,3] , ['B',3]) . coxeter _diagrain() 

Coxeter diagram of rank 8 

• bjmatrix(self ) 

Returns the exchange matrix of self 

sage: QuiverMutationTypeC ['A' ,5] ) .b_matrix() 



/ 1 
-10-10 

1 1 
0-10 




\ 1 



/ 



sage: QuiverMutationTypeC ['A' ,3] , ['B' 



' ,3] ) .bjnatrixO 
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• staiidard_quiver (self ) 

Returns the standard quiver of self. 

• cartanjmatrix(self ) 

Returns the Cartan matrix of self which is obtained from its exchange 
matrix by replacing the positive entries by negative, and replace the O's on 
the main diagonal by 2's. 

sage: QuiverMutationTypeCA' ,5) . cartan_matrix() 

/ 2 -1 \ 
-1 2-1 
0-1 2-1 
0-1 2-1 
V -1 2 / 
sage : QuiverMutationTypeC ['A' ,3] , ['B' ,3] ) . cartanjnatrixO 
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class_size (self ) 

Returns the number of quivers which are mutation-equivalent to self, up 
to isom orph ism (Warning: several class sizes are only conjectured, see 
Section |43|. 



sage : QuiverMutationType( ['A' , 22] , ['BD' , 16, 1] ) .class_size() 

4257164518523691840 
sage : QuiverMutationType( ['GR' , [4,9] ,3] ) . class_size() 

oo 

• dual (self) 

Returns the dual quiver mutation type of self. 

sage: QuiverMutationTypeCA' ,4) .dualO 

['A', 4] 

sage: QuiverMutationTypeCB' ,4) .dualO 

['C, 4] 

• is_irreducible (self ) 
Returns True, iff self is irreducible. 

sage: QuiverMutationTypeCA' ,4) . is_irreducible O 

True 

sage : QuiverMutationType( [' A' ,3] , ['B' ,3] ) . is_irreducible () 

False 

• is_mutation_f inite (self ) 

Returns True, iff self is of finite mutation type. 

sage : QuiverMutationType( ['GR' , [4,8] ,3] ) . is_mutation_f initeO 

True 

sage: QuiverMutationType( ['GR' , [4,9] ,3]) . isjnutation.f initeO 
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False 

• is_simply_laced(self ) 

Returns True, iff self is simply-laced. 

• is_f inite (self ) 

Returns True, iff self is of finite type. 

• is_af fine (self ) 

Returns True, iff self is of affine type. 

• is_elliptic (self ) 

Returns True, iff self is of elliptic type. 

• irreducible.components (self ) 

Returns a tuple containing the irreducible components of self. 

sage: QuiverMutationTypeCA' ,5) . irreducible_components () 

(['A', 4],) 

sage: QuiverMutationType( ['A' ,3] , ['B',3]) . irreducible_components () 
(['A', 3], ['B', 3]) 

• properties (self ) 

Prints all properties of self. See Section |4] for examples. 

6.3. Quiver. The next class we want to describe it the class Quiver. It allows 
numerous ways to construct a quiver, several examples were described in Section [3] 

• QuiverMutationType 

• list or tuple representing a quiver mutation type 

• ClusterSeed 

• matrix: a skew-symmetrizable matrix which represents the exchange ma- 
trix 

• Quiver 

• DiGraph: the digraph must represent a quiver 

• list of tuples representing the edge list of a digraph for a quiver 
The attributes of Quiver are given by 

• _M 

The exchange matrix of self. 

• jn 

The number of cluster variables (which is the number of columns in the 
exchange matrix). 

• _m 

The number of frozen variables (which is the number of rows minus the 
number of columns in the exchange matrix). 

• -description 

The string representation of self. 
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• _mutation_type 

The mutation type of self, if known, None otherwise. 

The methods of Quiver are given by 

• __init__(self , data, frozen=0) 

Frozen sets the later vertices to be frozen 

sage: Ql = Quiver ( [(0, 1) , (1 ,2) , (2,3)] ) ; Ql 

Quiver on 4 vertices 

sage: Q2 = Quiver([(0,l) , (1,2) , (2,3)] ,frozen=l) ; Q2 

Quiver on 4 vertices with 1 frozen vertex 

sage: Ql.showO 
sage: Q2.show() 




• __eq__(self, other) 

Returns True, iff the b-matrices of self and other coincide 

sage: Q = Quiver( ['A' ,5] ) 

sage: T = Q.mutate( 2, inplace=False ) 

sage: Q.__eq__( T ) 

False 

sage: T. mutate ( 2 ) 
sage: Q.__eq__( T ) 

True 

• _repr_(self) 

Returns the string representation of self 

sage: Q = Quiver( ['A' ,5] ) 
sage: Q._repr_() 

"Quiver on 5 vertices of type ['A', 5]" 

• plot (self , circular=False , directed=True , mark=None) 

Returns a random/circular and directed/undirected plot of self with a given 
vertex marked. 

• show(self, fig_size=l, circular=False , directed=True , mark=None) 

Shows the plot of self. 



GREGG MUSIKER AND CHRISTIAN STUMP 



interact (self , fig_size=l, circular=True) 

Starts an interactive mode, as shown in Figure [T] at the end of Section [3j 

save_image (self , filename, circular=False) 

Saves the plot of self to filename. 

bjnatrix(self ) 

Returns the exchange matrix of self 
sage: QuiverC [' A' ,4] ) .bjnatrixO 
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sage: Quiver ( ['B' ,4] ) .b_matrix() 
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sage: Quiver ( ['D' ,4] ) .b_matrix() 
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sage: Quiver(QuiverMutationType( [ [' A' ,2] , ['B' ,2]])) .bjnatrixO 
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digraph (self ) 

Returns the underlying digraph of self 
sage: QuiverC ['A' ,4] ) .digraphO 

Digraph on 4 vertices 

n(self) 

Returns the number of free vertices of self 

sage: Q = Quiver([(0,l) , (1,2) , (2,3)] ,frozen=l) 
sage: Q.nO 

3 

m(self ) 

Returns the number of frozen vertices of self 

sage: Q = Quiver([(0,l) , (1,2) , (2,3)] ,frozen=l) 
sage: Q.mO 

1 
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• cajionical-label (self , certif y=False) 

Returns an isomorphic quiver with canonical vertex labeling. This is based 
on the canonical labeling of digraphs using the corresponding method for 
digraphs by R.L. Miller based on [McKSlj . If certify is True, a dictionary 
of the relabeling is also returned 

sage: Quiver(['A' ,4]) . canonical_label (certif y=True) 

(Quiver on 4 vertices of type ['A', 4], {0:0,1:3,2:1,3:2}) 

• is_acyclic (self ) 

Returns True, iff self is acyclic. 

• is_bipartite (self , return_bipartition=False) 

Returns True, iff self is bipartite, if return_bipartition is True, the biparti- 
tion is returned 

sage: Quiver(['A' ,4]) . is_bipartite(return_bipartition=True) 
(set([0, 2]), set([l, 3])) 

• principal_restriction(self ) 

Returns the principal restricting of self. This is obtained from self by delet- 
ing all frozen variables. 

• principal_extension(self ) 

Returns the principal extension of self. This can be used only for seeds 
without frozen variables. Returns a new seed with exchange matrix of size 
2n X n given by the exchange matrix of self of size nx n with an additional 
identity matrix added below. 

• mutate (self , data, inplace=True) 

Mutates at a vertex or at a list of vertices, if inplace is True, self is modi- 
fied, otherwise a new quiver is returned. 

• mutation_sequence (self , sequence, show_sequence=False , 

f ig_size=l . 2) 

Returns a list of quivers obtained from a sequence of mutations. If the pa- 
rameter show_sequence is True, the sequence is shown with a given fig_size. 

• reorient (self , data) 

Reorients self with respect to the given total order, or with respect to an 
iterator of edges in self to be reverted. 

Warning: This often will change the mutation class of the quiver except 



if the quiver is a tree (see Theorem 4.7). 



mutation_class_iter (self , depth=inf inity , 
show_depth=False , return_paths=False , 

data_type=' quiver' , up_to_equivalence=True , 
only_sink_source=False) 

Returns an iterator which goes through the mutation class of self depending 
on several parameters 
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— depth: integer, only quivers with distance at most depth from self are 
returned 

— show_depth: if True, the actual depth of the mutation is shown 

— return_paths: if True, a shortest path of mutation sequences from 
self to the given quiver is returned as well 

— data_type; can be one of the following: 

quiver, matrix, digraph, dig6, path 

— up_to_equivaleiice: if True, only quivers up to equivalence are con- 
sidered 

— sink_source: if True, only mutations at sinks and sources are applied 

mutatioii_class (self , depth=inf inity , 
show_depth=False , return_paths=False , 

data_type=' quiver' , up_to_equivalence=True , 
only_sink_source=False) 

Returns a list of all quivers in the corresponding iterator. 

group_of jnutations (self ) 

Returns the group of mutations of self. Warning: The permutation group 
is very big! This group differs for quivers and for cluster seeds, as different 
cluster seeds may have the same exchange matrix and thus the same quiver. 
This group is defined to be the group of permutation given as follows. The 
ground set is the mutation class of self without taking equivalence of quivers 
into account, and the group is generated by the n involutions on this set 
given by mutation at the n different vertices. Observe that the analogous 
operation on the mutation class up to equivalence does not give a group 
(this can be easily checked in type A^). Basically nothing is known about 
this group. 

sage: Q = Quiver( ['A' ,2] ) 
sage: Q.group_of _mutations() 

Permutation Group with generators [(1,2)] 

sage: Q = Quiver( [' A' ,3] ) 
sage: Q. group _of_mutations() 

Permutation Group with generators 
[(1,2) (3, 4) (5, 9) (6, 7) (8, 12) (10, 11) (13, 14), 

(1.3) (2,5) (4,6) (7,14) (8,11) (9,13) (10,12) , 

(1.4) (2,3) (5, 10) (6,8) (7,13) (9,14) (11,12)] 

sage: Q = Quiver( ['B' ,2] ) 
sage: Q.group_of .mutations () 

Permutation Group with generators [(1,2)] 

sage: Q = Quiver( ['B' ,3] ) 
sage: Q.group_of_mutations() 

Permutation Group with generators [(1 ,2) (3 ,4) (5 ,6) (7 , 10) (8 ,9) , 
(1,3) (2, 6) (4, 5) (7, 9) (8, 10), (1 ,4) (2 ,3) (5 ,7) (6, 8) (9 , 10)] 

sage: Q = Quiver( ['A' , 1] ) 

sage: Q.group_of_mutations() . cardinality () 

1 



sage 



Q = Quiver( ['A' ,2] ) 
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sage: Q.group_of .mutations () . cardinality () 

2 

sage: Q = QuiverC ['A' ,3] ) 

sage: Q.group_of jnutationsO . cardinality () 

322560 

• is_f inite (self ) 

Returns True, iff self is of finite type. This is done by checking if it is 
mutation-equivalent to a quiver of finite type. 

• is_tiiutation_f inite (self , nr_of _checks=None , return_path=False) 
Returns True, iff self if of finite mutation type. Warning: The algorithm is 
non-deterministic and uses random mutations in various directions. Might 
theoretically result in a wrong True return. The number of checks can be 
set, the default is 1000 times the number of vertices of self. If return_path 
is True, then a path to a non- mutation- finite quiver is returned, if found. 

• mutation_type(self ) 

Returns the mutation type of self if it can be determined. 

— First, it is checked if self is mutation-equivalent to a quiver of a classical 
type using the descriptions of the classification types, 

— then, it is checked if self is contained in an exceptional mutation class 
which are hard-coded, 

— if it was not possible to determine the mutation type, it is checked if 
self is mutation-finite or infinite 

Warning: The algorithm to determine quivers of mutation type £>„ (which 
is ['D' ,n,l]) is not yet implemented! 

6.4. ClusterSeed. The constructor of the class ClusterSeed allows the same in- 
put as the class Quiver to construct a cluster seed. Moreover, many attributes and 
methods for cluster seeds and for quivers coincide. Often, the cluster seed simply 
calls the quiver method. 

• QuiverMutationType 

• list or tuple representing a quiver mutation type 

• ClusterSeed 

• matrix: a skew-symmetrizable matrix which represents the exchange ma- 
trix 

• Quiver 

• DiGraph: the digraph must represent a quiver 

• list of tuples representing the edge list of a digraph for a quiver 
The attributes of ClusterSeed are given by 

• J4 

The exchange matrix of self. 

• .cluster 

The cluster as a list of cluster variables. 

• jn 

The number of cluster variables (which is the number of columns in the 
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exchange matrix). 



• _m 



The number of frozen variables (which is the number of rows — the number 
of columns in the exchange matrix). 

• _R 

The base ring in which the cluster variables live. 

• _quiver 

The quiver attached to self. 

• .description 

The string representation of self. 



• _mutation_type 

The mutation type of self, if known, None otherwise 

The methods of ClusterSeed are given by 

• __init__(self , data, frozen=0) 

Frozen sets the later vertices to be frozen 

sage: SI = ClusterSeedC [(0,1) , (1,2) , (2,3)] ) ; SI 

A seed for a cluster algebra of rank 4 
sage: Q2 = Quiver ( [(0, 1) , (1 ,2) , (2,3)] ,frozen=l) ; Q2 

A seed for a cluster algebra of rank 3 with 1 frozen variable 
sage: Ql .bjnatrixO ; Q2 .b_matrix() 



/O 1 00\ /O 1 0\ 

-10 10 -10 1 

0-101 0-10 

\0 0-10/ \0 0-1/ 



• __eq__(self, other) 

Returns True, iff self and other have the same exchange matrix and the 
same cluster. 



• _repr_(self) 

Returns the string representation of self 

sage: S = ClusterSeed( [' A' ,3] ) ; S._repr_() 

"A seed for a cluster algebra of rank 3 of type ['A', 3]" 

• plot (self , circular=False , mark=None) 

Returns a random/circular plot of self with a given marked vertex. Calls 
the method for quivers. 



• show(self, fig_size=l, circular=False , mark=None) 

Shows the plot of self. 



A COMPENDIUM ON THE CLUSTER ALGEBRA AND QUIVER PACKAGE IN Sage 



• interact (self , fig_size=l, circular=True) 

Starts an interactive mode, as shown in Figure [T] at the end of Section 

• save_image (self , filename, circular=False) 

Saves a plot of self to filename. 

• bjnatrix(self ) 

Returns the exchange matrix of self. 

• cluster (self ) 

Returns the cluster of self 

sage: S = ClusterSeed( ['A' ,3] ) ; S. cluster () 

[xo, 3:1,3:2] 
sage: S.inutate(O) ; S.clusterO 

'xi + 1 

,Xi,X2 

Xo 

sage: S. mutate ( 1) ; S.clusterO 



xi + 1 X0X2 + a;i + 1 

, ,X2 

Xo XoXi 

• ground_f ield(self ) 

Returns the ground field in which the cluster variables of self live 

sage: S.ground_field() 

Fraction Field of Multivariate Polynomial Ring 
in xO, xl, x2 over Rational Field 

• x(self,k) 

Returns the fcth initial cluster variable of self. 

• y(self ,k) 

Returns the fcth frozen variable of self. 

• n(self) 

Returns the number of cluster variables of self. 

• m(self) 

Returns the number of frozen variables of self. 

• exchangeable.variables (self ) 

Returns a list of all cluster variables of self. 

• f rozen_variables (self ) 

Returns a list of all frozen variables of self. 
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quiver (self ) 

Returns the Quiver associated to self. 

is_acyclic (self ) 

Returns True, iff self is acyclic. 

is_bipartite (self , return_bipartition=False) 

Returns True, iff self is bipartite, if return_bipartition is True, the biparti- 
tion is returned 

sage: ClusterSeedC [' A' ,4]) . is_bipartite (return_bipartition=True) 
(set([0, 2]), set([l, 3])) 

mutate (self , sequence, inplace=True) 

Mutates at an index or at a list of indices, if inplace is True, self is modified, 
otherwise a new cluster seed is returned. 

mutation_sequence (self , sequence, show_sequence=False , 

f ig_size=l . 2,return_output=' seed' ) 
Returns a list depending on return_output obtained from a sequence of 
mutations. If show -sequence is True, the sequence is shown with a given 
fig_size. The possible outputs are 

— 'seed': a list of cluster seeds is returned 

— 'matrix' : a list of exchange matrices is returned 

— ' var ' : a list of cluster variables is returned 

principal_restrict ion (self ) 

Returns the principal restriction of self. This is obtained from self by delet- 
ing all frozen variables. 

principal_extension(self ) 

Returns the principal extension of self. This can be used only for seeds 
without frozen variables. Returns a new seed with exchange matrix of size 
2n X n given by the exchange matrix of self of size nx n with an additional 
identity matrix added below. 

reorient (self , data) 

Reorients self by reorienting the corresponding quiver. Calls the method 
for quivers. 

set_cluster (self , cluster) 

Sets the set of clusters of self to cluster. 

reset .cluster (self ) 

Sets the set of clusters of self back to the initial cluster. 

mutation_class_iter (self , depth=inf inity , 
show_depth=False , return_paths=False , 

up_to_equivalence=True , only_sink_source=False) 
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Returns an iterator which goes through the mutation class of self depending 
on several parameters 

— depth: integer, only quivers with distance at most depth from self are 
returned 

— show_depth: if True, the actual depth of the mutation is shown 

— return_paths: if True, a shortest path of mutation sequences from 
self to the given quiver is returned as well 

— up_to_equivalence: if True, only quivers up to equivalence are con- 
sidered 

— only_sink_source: if True, only mutations at sinks and sources are 
applied 

• mutation_class (self , depth=inf inity , 

show_depth=False , return_paths=False , 

up_to_equivalence=True , only_sink_source=False) 

Returns a list of all quivers in the corresponding iterator 

• cluster_class_iter (self , depth=inf inity , show_depth=False , 

up_to_equivalence=True) 

Returns an iterator through all clusters mutation-equivalent to self up to 
a given depth. Moreover, it is possible to show the actual depth together 
with several parameters, or to output clusters as labeled seeds. 

• cluster_class(self , depth=inf inity , show_deptli=False , 

up_to_equivalence=True) 

Returns a list of all clusters mutation-equivalent to self up to a given depth. 
Moreover, it is possible to show the actual depth together with several pa- 
rameters, or to output cluster as labeled seeds. 

• bjiiatrix_class_iter (self , deptli=inf inity , 

up_to_equivalence=True) 

Returns an iterator through all matrices mutation- equivalent to self up to 
a given depth, and up to permutation of rows and columns unless specified 
otherwise. 

• bjnatrix_class (self , depth=inf inity , up_to_equivalence=True) 

Returns a list of all matrices mutation-equivalent to self up to a given depth, 
and up to permutation of rows and columns unless specified otherwise. 

• variable_class_iter (self , depth=inf inity , 

ignore_bipartite_belt=False) 

Returns an iterator through all variables obtained from self by mutations 
up to a given depth. Warning: If at some point a bipartite seed is reached, 
another algorithm is used unless the parameter ignore_bipartite_belt is 
set to be True. See the description in Section [44} 

• variable_class (self , depth=inf inity , 

ignore_bipartite_belt=False) 
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Returns a list of all variables obtained from self by mutations up to a given 
depth. Warning: If at some point a bipartite seed is reached, another 
algorithm is used unless the parameter ignore_bipartite_belt is set to 



be True. See the description in Section 4.4 



group_of jnutations (self ) 

Returns the group of mutations of self. Warning: The permutation group 
is very big! This group differs for quivers and for cluster seeds, as different 
cluster seeds may have the same exchange matrix and thus the same quiver. 
This group is defined to be the group of permutation given as follows. The 
ground set is the mutation class of self without taking equivalence of seeds 
into account, and the group is generated by the n involutions on this set 
given by mutation at the n different vertices. Observe that the analogous 
operation on the mutation class up to equivalence does not give a group 
(this can be easily checked in type A^). Basically nothing is known about 
this group. 

sage: S = ClusterSeedC ['A' ,2] ) 
sage: S.group_of_mutations() 

Permutation Group with generators [(1 ,2) (3 ,4) (5 ,6) (7 ,9) (8 , 10) , 
(1,3)(2,5)(4,7)(6,8)(9,10)] 

sage: S = ClusterSeedC ['B' ,2] ) 
sage: S.group_of_mutations() 

Permutation Group with generators 
[(1,2)(3,4)(5,6), (1,3)(2,5)(4,6)] 

sage: Q = ClusterSeedC ['A' , 1] ) 

sage: Q.group_of_mutations() . cardinality () 



sage: Q = ClusterSeedC ['A' ,2] ) 

sage: Q.group_of_mutationsC) . cardinality C) 

10 

sage: Q = ClusterSeedC ['A' ,3] ) 

sage: Q.group_of .mutations C) . cardinality C) 

705438720 



is_f inite (self ) 

Returns True, iff self is of finite type. Calls the method for the quiver of self. 

isjnutation_f inite (self , nr_of _checks=None , return_path=False) 
Returns True, iff self is of finite mutation type. Calls the method for the 
quiver of self. 

mutation_type(self ) 

Returns the mutation type of self, if possible. Calls the method for the 
quiver of self. 
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6.5. ClusterVariable. By definition, a cluster variable is an element in the field 
of rational function in n variablef]^ The class ClusterVariable provides two extra 
features for cluster variables: 

(1) The connection to almost positive roots in finite types (positive roots are 
not yet provided in sage for affine types). 

(2) An ordering for cluster variables which is inspired by its connection to 
almost positive roots: 

• They are ordered first by total degree of the denominator (in particu- 
lar, the variables in the initial seed come first in natural order), 

• If the degree is equal and positive, they are ordered lexicographically 
with a;o > Xl > . . . > Xn-i- 

sage: for f in ClusterSeedC ['A' ,2] ) . variable_class() : 
. . . . : print f , f . almost_positive_root () 



xo — ai 

Xl ~ a2 

{xi + l)/xo ai 

(a::o + 02 

(xo + a;i + l)/(a::oa;i) ai + 02 



Two further examples of the ordering can be found in Section 4.2 It is planned to 



include more functionalities for the cluster variable class in the future. 
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